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just  like  shocks  created  by  more  traditional  methods.  In  one  experiment  we  examined  laser-produced-shock  solid-surface 
interactions  and  observed  expected  phenomena  such  as  Marh  stems  and  triple  points,  and  also  measured  a  new  phenomenon 
termed  a  blast  wave  decursor.  In  second  experiment  we  found  that  shocks  become  unstable  if  they  propagate  through  a  gas 
which  has  a  low  adiabatic  index  and  we  measured  the  growth  rate  of  die  instability.  In  a  third  experiment  we  have  shown  that 
a  high  Mach  number  shock  dramatically  enhances  the  structure  of  a  turbulent  field  through  which  it  passes  and  that  the  shock 
is  itself  badly  distorted.  This  result  is  unexpected  since  common  wisdom  has  it  that  high  Mach  number  shocks  would  self-heal 
as  they  pass  through  a  turbulent  field. 
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EXPERIMENTAL  STUDIES  OF  VERY-HIGH 
MACH  NUMBER  HYDRODYNAMICS 


J.  Grun*,  C.  K.  Manka1,  B.  H.  Ripin1,  A-  C.  Buckingham2,  and  I.  Kohl  berg3 
1  Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington  DC  20375,  USA 
^Lawrence  Livermore  National  Laboratory,  Livermore,  Ca,  USA 
3Kohlberg  Associates,  Inc.,  Alexandria,  VA  22304,  USA 

Abstract.  We  present  results  of  experiments  on  very-high  Mach  number  (>100)  shocks  and  veiy-high 
Mach  number  and  Reynolds  number  (>100, 10^)  turbulence.  Such  high  Mach  number  hydrodynamics  are 
initiated  with  a  powerful  laser  pulse  driver.  We  show  that  shocks  created  with  a  laser  driver  follow  the 
Taylor-Sedov  self-similar  solution  and  scale  via  the  Sachs  scaling  law  just  like  shocks  created  by  mote 
traditional  methods.  In  one  experiment  we  examined  laser-produced-shock  solid-surface  interactions  and 
observed  expected  phenomena  such  as  Mach  stems  and  triple  points,  and  also  measured  a  new 
phenomenon  termed  a  blast  wave  decursor.  In  second  experiment  we  found  that  shocks  become  unstable 
if  they  propagate  through  a  gas  which  has  a  low  adiabatic  index  and  we  measured  the  growth  rate  of  the 
instability.  In  a  third  experiment  we  have  shown  that  a  high  Mach  number  shock  dramatically  enhances 
the  structure  of  a  turbulent  field  through  which  it  passes  and  that  the  shock  is  itself  is  badly  distorted. 
This  result  is  unexpected  since  common  wisdom  has  it  that  high  Mach  number  shocks  would  self-heal  as 
they  pass  through  a  turbulent  field. 


Key  words:  Laser  matter  interaction.  Shock,  Turbulence,  Blast  wave.  Decursor,  Shock  instability. 
Shock  turbulence  interaction 


1.  Introduction 

High  Mach  number  and  high  Reynolds  number  hydrodynamics  occurs  in  man  made  and  natural 
phenomena.  Examples  include  wakes  and  shocks  generated  by  hypersonic  vehicles,  shocks  from  surface 
discontinuities  on  hypersonic  vehicles  interacting  with  a  turbulent  boundary  layer  near  die  vehicle 
surface,  shocks  induced  by  jet  or  slot  coolant  injection  interacting  with  separated  shear  layers,  flows 
within  ramjet  engines,  and  supemovae  explosions.  The  nature  of  the  interactions  between  the  different 
flows  can  influence  eddy  transport,  component  mixing,  diffusion,  surface  shear,  and  heat  transfer  •  and 
thus  have  important  consequences  in  practical  systems.  Shock  instabilities  are  thought  to  play  a  role  in 
structures  seen  in  supemovae  remnants  and  in  galaxy  and  star  formation. 

Very  high  Mach  number  and  Reynolds  number  flows  are  initiated  at  the  Naval  Research  Laboratory 
(NRL)  using  a  powerful  laser  pulse  driver.  In  these  experiments  a  millimeter-diameter,  few-micron 
thick  target  is  placed  in  an  evacuated  chamber,  which  is  then  backfilled  with  gas.  The  target  is  irradiated 
and  super-heated  at  an  intensity  of  1  to  100  Terawatts/cm2  by  a  focused  pulse  from  the  Pharos  in  Nd- 
glass  laser,  which  has  three  beams  and  produces  up  to  1500  Joules  per  pulse.  A  very  thin  (-  1  micron) 
layer  ablates  from  the  target  surface,  expands  rapidly  into  the  background  gas,  and  -  much  like  the 
products  of  a  chemical  explosion  -  forms  a  shock  wave  expanding  in  the  direction  of  the  incident  laser 
pulse  (Gnm,  et  al  - 1981,1986,  Ripin,  et  al  - 1986,  Stamper,  et  al  - 1988).  The  rest  of  the  target  is  pushed 
into  the  background  gas  by  the  rocket-like  effect  of  the  ablating  material.  This  part  of  the  target  rapidly 
falls  apart  due  to  the  Raleigh  Taylor  instability  (Gnm,  et  al  -  1984,  1987)  creating  high  Mach  number 
and  high  Reynolds  number  turbulence  in  the  ambient  gas. 

Shock  waves  created  by  a  laser  pulse  can  have  Mach  numbers  up  to  a  few  hundred,  pressures  up  to  tens  of 
megabar,  and  temperatures  from  a  fraction  of  eV  to  greater  than  1  KeV.  Turbulence  spatial  structure  can 
have  Reynolds  numbers  of  the  order  of  one  to  one  hundred  million.  By  using  high  atomic  number  target 
material  x-ray  radiation  of  a  few  KeV  can  be  produced  together  with  the  hydrodynamics.  Typical  spatial 
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Figure  1.  Setup  of  a  typical  experiment:  Ablation  from  a  target  surface  drives  a  high  Mach  number 
shock  into  a  background  gas  while  the  rest  of  the  target  generates  turbulence.  The  shock  front  or 
turbulence  structures  are  recorded  with  optical  shadowgraphy  and  interferometry  techniques  that  are 
sensitive  to  index-of-iefraction  fluctuations.  Examples  of  such  methods  are  bright-fteld  shadowgraphy, 
dark  field  shadowgraphy,  phase-coo trast-shadowgraphy,  shearing  interferometry  and  holographic 
interferometry.  In  our  experiments  up  to  four  sequential  frames  can  be  recorded  with  these  methods.  The 
last  three  methods  produce  quantitative  results  that  can  be  reliably  compared  to  theory.  In  addition,  the 
visible  emission  from  the  blast  front  is  photographed  with  two  very  fast  (100-psec  shutter),  four-frame, 
microchannel-plate  intensifier  cameras.  Thus,  on  each  shot  we  can  record  four  shadowgraph  or 
interferometric  images,  and  eight  emission  frames,  getting  twelve  photographs  of  the  blast  wave  at 
different  times  in  its  evolution.  Gas  emission  spectra  are  measured  with  time  and  space  resolved 
spectroscopy.  Rom  these  measurements  the  state  of  the  gas  is  inferred.  Fiber-optic  sensors  can  be  used  to 
measure  shock  pressure. 
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and  temporal  scales  arecentimeter  and  microsecond.  A  special  feature  of  the  laser  driver  technique  is  that 
the  ratio  of  energy  to  mass  is  an  extremely  high  - 109  joule/gram,  which  is  difficult  to  obtain  by  other 
means. 


In  this  paper  we  summarize  the  results  of  three  experiments:  (1)  an  experiment  which  shows  that  laser- 
produced  shocks  behave  like  shocks  produced  by  other  methods.  We  demonstrate  the  formation  of  Mach 
stems  and  triple  points,  and  show  a  new  phenomenon  termed  a  blast  wave  decursor.  (2)  an  experiment 
which  demonstrates  that  shocks  become  unstable  if  they  propagate  through  a  gas  which  has  a  low 
adiabatic  index,  and  (3)  an  experiment  which  shows  that  a  high  Mach  number  shock  enhances  the 
structure  of  a  turbulent  field  through  which  it  passes  and  is  itself  distorted. 


2.  Experiments 

2.1  Laser-produced  shocks  behave  normally 

We  proceed  to  show  that  laser-produced  shocks  obey  the  Taylor  -  von  Neumann  -  Sedov  relation 
(Zel'dovich  and  Raizer  - 1966)  and  the  Sachs  scaling  law  (Sachs  •  1944).  Ideal  shocks  (blast  waves) 
propagate  according  to  the  relation 

R  =  C(E/p0)l/s  t275.  wilh  c  =  (1) 

where  E  is  the  explosion  energy,  p0  is  the  ambient-gas  density,  t  is  the  time  of  observation,  and  y  is  the 
adiabatic  index  or  ratio  of  specific  heats.  These  expressions  assume  spherical  symmetry,  energy  and 
momentum  conservation,  and  an  instantaneous,  mass-less,  point  explosion.  In  our  experiments  we  varied 
the  laser  energy,  laser  focal-spot  size  (and  therefore  the  ablation  velocity),  pressure,  and  observation 
times.  We  found  that  as  long  as  the  ambient  pressure  is  high  enough  (>0.5  ton)  to  place  the  experiment 
in  a  collisional  regime,  laser-produced  shocks  follow  equation  (1).  From  the  data  and  equation  (1)  we 
determined  that  £  =  1.0  ±  0.1  in  a  nitrogen  ambient  gas,  which  implies  a  y  of  1.3  (Ripin  et  al  1986). 
These  results  are  shown  in  Figure  2. 

Shock  overpressure  is  determined  from  dark-field  photographs  such  as  that  in  Fig.  2  using  the  expressions 
(Kinney  et  al  -  1985): 


(dR/dt) 

"  c0 


(3) 


=  [3(y- 1  )(Y+-l)2/4K(3y- 1)] 1/2  [E/pJ  R*2 


(4) 


where  P  is  the  peak  over-pressure,  P0  is  the  ambient  pressure,  M  is  the  Mach  number,  and  Cq  is  the  sound 
speed  in  the  ambient  gas.  Figure  3  plots  the  overpressures  inferred  from  equations  2-4  as  a  function  of 
distance  from  the  explosion  scaled  using  the  Sachs  scaling  law.  This  law,  which  is  based  on  geometrical 
similarity  principles  and  conservation  of  momentum,  is  used  to  compare  experiments  performed  with 
different  explosive  yields  and  at  different  ambient  pressures.  Also  plotted  are  curves  expected  of  a 
Taylor  •  voo  Neumann  •  Sedov  shock  and  shocks  created  with  explosives  (Kinney  et  al  -  1985).  Figure  3 
shows  that  collisional  laser-produced  shocks  follow  the  ideal  scaling  law  and  continue  to  do  so  even  where 
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Figure  2.  Top:  Typical  dark-field  shadowgraph  of  a  shock  produced  by  irradiating  and  aluminum  target. 
Bottom:  A  graph  of  shock  location  versus  time  shows  that  laser-produced  shocks  obey  the  Taylor  -  von 
Neumann  -  Sedov  relation . 
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Figure  3.  Overpressures  of  laser-created  shocks  compared  to  Taylor  -  von  Neumann  -  Sedov  shocks  and 
shocks  created  by  chemical  explosions.  Data  is  compiled  on  one  plot  utilizing  the  Sachs  scaling  law 
(Gran,  et  aL  1991). 
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Figure  4.  Decursor  shock  near  the  surface  of  an  aluminum  plane.  Note  the  triple  points  and  Mach  stem 
in  the  reflected  shock. 
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shocks  produced  with  chemical  explosives  deviate  from  ideal  scaling.  Furthermore,  the  laser  experiment 
produces  shocks  with  overpressures  greater  than  100  kbar,  which  is  more  than  1000  times  greater  than  the 
capability  of  conventional  explosives.  All  this  is.  a  consequence  of  the  very  large  energy  to  mass  ratio  of 
laser-produced  "explosions". 

It  is  well  known  that  a  blast  wave  colliding  with  a  planar  surface  sets  up  a  reflection  that  propagates  back 
into  the  ambient  gas  heated  by  the  passage  of  the  incident  front  The  interacting  incident  and  reflected 
blast  fronts  form  Mach  stems  and  triple  points.  If  the  sound  speed  near  the  reflecting  surface  is  higher 
than  in  the  rest  of  die  ambient  gas,  part  of  the  blast  front  moves  faster  and  outruns  the  bulk  of  the  blast 
wave  forming  a  structure  called  a  shock  precursor  (Glasstooe  and  Dolan  •  1977,  Giowacki,  et  al  -1986). 
When  the  sound  speed  near  the  reflecting  surface  is  slower  than  in  the  rest  of  the  ambient  gas  a  part  of  the 
blast  wave  can  lag  behind  the  rest  The  existence  of  such  lagging  structure,  called  a  decursor,  has  been 
predicted  theoretically  and  was  observed  first  in  our  experiments.  Figure  4  shows  examples  of  Mach 
stems,  triple  points,  and  decursors  measured  in  our  experiments.  See  Grun,  et  al.  1991  for  a  discussion  of 
these  phenomena. 


22  High  Mach  number  shocks  become  unstable  in  a  low  y  gas 

Shock  instabilities  are  thought  to  contribute  to  die  structuring  Observed  in  supernovae  and  play  a  role  in 
the  formation  of  stars  and  pi”"*  (Chevalier  -  1976),  Gerola  and  Seiden,  -1978,  Ostriker  and  Cowie  - 
1981),  Trimble  -  1988).  Theories  of  shock  instability  in  uniform  ambient  have  been  accompanied  by 
considerable  controversy  (Isenberg  -  1977,  Cheng  -  1979,  Bernstein  and  Book  -  1980,  Gaffet-  1984, 
Kohlberg  -1989,  Newman  -  1980,  Vishniac  -  1983),  Ryu  and  Vishniac  -  1987,  Vishniac  and  Ryu  -1989). 
with  some  concluding  that  such  instabilities  do  not  exist  (Newman  -  1980).  Laser  produced  shocks  have 
characteristics  similiar  in  some  respects  to  supernovae  ( Ripin,  et  al  - 1990).  Our  experiments  (Grun  et  al  - 
1991)  show  conclusively  that  shocks  are  unstable  if  they  propagate  through  a  uniform  gas  with  a  low 
adiabatic  index  y  (such  as  *f=1.06).  These  results  were  consistent  with  the  1989  theory  of  Vishniac  and 
Ryu.  A  nonlinear  analysis  of  our  results  was  performed  by  Low  and  Norman  in  1992. 

The  experiments  were  designed  to  compare  the  propagation  of  shocks  through  two  gases:  nitrogen  and 
xenon.  A  basic  difference  between  these  gases  is  that  nitrogen  has  an  adiabatic  index  Yn=13  ±0. 1,  while 
xenon  has  a  much  lower  adiabatic  index  ^=1.06  ±  0.02.  To  determine  ife  we  utilized  measurements 
which  showed  that  shocks  in  both  nitrogen  and  xenon  propagate  according  to  equation  (1).  Therefore,  by 
dividing  the  measured  radius  R  in  nitrogen  by  the  measured  radius  R  in  xenon  at  any  given  time  we 
arrived  at  a  relationship  between  y^  yN,  and  foe  mass  of  each  gas  species.  Solving  for  y^  as  a  function 
of  yN,  we  found  that  as  -ft,  varies  from  1  to  5/3,  y^  varies  from  1  to  1.13.  Hence,  for  any  reasonable  value 
of  yN  foe  value  of  y^  must  be  less  than  1.13.  Moreover,  since  we  have  determined  Yn  to  **  *-3  -0.1, 
Yxc  must  be  1.06  ±  0.02. 

Y&  is  lower  than  Yn  because  xenon  gas  radiates  much  more  than  nitrogen  gas.  Radiation  increases  foe 
degrees  of  freedom  within  a  gas  and  hence  reduces  its  effective  y.  This  was  demonstrated  by  examining 
foe  spectrum  of  light  emanating  from  a  point  some  distance  in  front  of  foe  laser’s  focal  spot  In  nitrogen 
gas  foe  laser-induced  explosion  produces  little  measurable  emission  prior  to  foe  arrival  of  foe  shock. 
Immediately  before  foe  shock  arrives  at  foe  observation  point  there  is  a  slight  increase  of  N1+  and  N2+ 
lines:  These  lines  are  probably  excited  by  UV  or  beat  from  foe  shock.  When  the  shock  front  reaches  foe 
observation  point  there  is  a  sudden  increase  in  emission  from  foe  N1+,  N2+,  and  target  C2*  lines,  as  well 
as  an  inrr<»a<a»  in  continuum  emission.  In  contrast  xenon  radiates  copiously  in  many  XeI+  and  some 
Xe2+  lines  from  foe  moment  foe  laser  bits  foe  foil  and  well  before  foe  shock  reaches  foe  observation  point 
We  conclude,  therefore,  that  it  is  foe  radiation  in  xenon  which  reduces  its  effective  y  below  that  of 
nitrogen. 


7 


time  (ns)  log(2*FVL) 


Figure  5.  Top:  Dark-field  shadowgraph  of:  a  stable  blast  wave  L.  nitrogen  gas,  and  an  unstable  blast 
wave  in  xenon  gas.  Bottom:  Instability  growth  amplitude  as  a  function  of  time  for  different  values  of 
log(kR)and  Growth  exponent  S(kR). 


8 


Imaging  diagnostics  show  that  shocks  propagating  through  nitrogen  gas  are  stable  and  smooth.  In 
startling  contrast,  shocks  propagating  through  xenon  gas  become  wrinkled  like  a  dried  prune!  This 
wrinkling  was  quantified  by  tracing  and  then  Fourier  transforming  the  outer  edge  of  the  shock  front  and 
presenting  the  result  as  A,(k)/R  versus  k>g(kR),  where  A,(k)  is  the  full  amplitude  of  the  mode  with 
wavenumber  k  at  time  t,  and  R  is  the  average  radius  of  the  shock  boundary.  A  power  law  of  the  form 

AJ(k)/R  s  was  fit  to  the  A,(k)/R  vs.  time  data  and  the  exponent  S(kR)  was  determined.  We  found 
growth  occurring  for  modes  satisfying  0.7  <log(kR)  <  2.  Maximum  growth  occurred  at  ?og(kR)  *  1  where 
S  *  1.6,  and  minimum  growth,  with  S  «  0.3,  occurred  at  logflcR)  =  2. 

2,3  High  Mach  number  shocks  amplify  turbulence  spectrum 

Consider  what  happens  when  a  shock  traverses  a  turbulent  flow.  To  zero'th  order  the  turbulence  can  be 
viewed  as  a  collection  of  positive  and  negative  "shock  lenses"  of  various  powers,  sizes,  and  dimensions.  A 
shock  passing  through  those  "lenses"  breaks-up  into  sbocklets  that  focus  in  scone  places  and  defocus  in 
others.  The  focused  sbocklets  quickly  become  nonlinear,  deposit  their  energy  in  the  turbulent  field,  and 
thereby  alter  its  structure.  Such  shock-turbulence  interactions  have  been  studied  both  computationally  and 
experimentally  for  low  Mach  number  shocks  (Trolier  and  Duffy  -  1985,  Hartung  and  Duffy  -  1986;  Smits 
and  Muck  -  1987;  Keller  and  Merzkirdb  -  1990;  Buckingham  -  1987,  1989,  1990,  1991;  Growley  and 
Burk  -  1991).  Significant  enhancement  of  the  turbulence  structure  have  been  observed.  However,  there  is 
a  commonly  held  view  that  in  the  hypersonic  regime  (Mach  >6)  reinforcement  of  shock  strength  by 
nonlinear  pressure  field  interactions  heals  any  shock  front  irregularities  and  suppresses  shock  distortion 
as  well  as  its  influence  on  the  turbulent  field. 

Our  current  work  involves  measuring  the  interaction  of  very  high  Mach  number  shocks  with  turbulent 
flows  (Gran,  et  al  -  1992,  Buckingham  and  Grun  •  1993).  Such  work  is  needed  because  information  on 
the  strongly  interactive  influence  of  high  Mach  number  (>6)  and  high  Reynolds  number  (>10^  )  is 
virtually  non  existent.  Experimental  data  at  Mach<6  in  high  Reynolds  Number  flows  are  neither  useful 
nor  extrapolateable  to  higher  flow  Mach  number  conditions  because  low  Mach  number  flows  do  not 
exhibit  sufficiently  strong  acoustical -density-temperature-vcloaty  coupling  to  allow  separation  and 
analysis  of  these  important  physical  effects.  Furthermore,  without  the  ability  to  vary  important 
parameters,  such  as  Mach  number,  over  a  broad  hypersonic  range  understanding  of  Mach-Reynolds 
number  coupling  is  difficult  The  NRL  experiment  vary  the  Mach  number  over  a  broad  hypersonic  range 
and  through  systematic  variation  of  this  parameter  we  will  examine  the  coupling  of  Reynolds  number 
and  Mach  number. 

Initial  results  indicate  that  contrary  to  common  expectations  a  high  Mach  number  (M-100)  shock  does 
enhance  the  structure  of  a  turbulent  field  through  which  it  passes.  In  addition,  the  shock  itself  is  badly 
distorted  as  it  passes  through  the  turbulence.  In  the  experiment  turbulent  flow  is  produced  by  accelerating 
a  thin  foil  into  a  5-Torr  N2  ambient  background  gas  by  irradiating  the  foil  with  one  300  J  pulse  from  the 
PHAROS  III  laser.  A  second  Pharos  in  beam  creates  a  high  Mach  number  shock  which  flows  over  the 
turbulent  flow  created  by  the  first  beam.  The  power  spectral  density  (PSD)  of  the  turbulent  flow,  which 
provides  a  quantitative  measure  of  the  turbulence,  is  measured  with  phase-contrast  microscopy.  Figure  6 
shows  the  experimental  result  and  a  comparison  of  experimental  pre-shocked  and  post-shocked  turbulence 
spectrum  with  predictions  by  Lawrence  Livermore  Laboratory's  large  eddy  simulation  code  (Buckingham 
and  Grun  -  1993). 

3.  Summary 

We  have  shown  that  a  kilojoule  laser  is  a  suitable  driver  for  launching  and  studying  very  high  Mach 
number  hydrodynamic  flows  •  flows  that  are  difficult  to  study  by  any  other  method.  Our  experiments 
produced  the  first  observations  of  a  shock  decursor  as  well  as  the  first  measurements  of  shock 
instability^nd  unexpected  measurements  of  the  intensification  of  a  turbulent  structure  by  a  very  high 
Mach  number  shock. 
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Figure  6:  Top:  A  high  Mach  number  shock  changes  the  structure  of  turbulence  through  which  it  passes. 
In  this  experiment  one  beam  of  the  Pharos  laser  created  a  structured  flow  while  another  created  a  high 
Mach  number  shock.  Laser-beam  and  target  geometries  were  arranged  so  that  the  shock  would  flow  over 
the  structure.  Dark-field  photographs  like  the  one  here  produced  a  striking  visual  record  of  the 
enhancement  in  the  pre-existing  structure  as  it  is  over-run  by  a  shock.  Bottom:  Comparison  of  pre  and 
post  shocked  turbulence  structure.  (The  experiment  and  code  differ  at  high  wave  number  because  the 
experiment  is  in  spherical  and  the  code  in  cylindrical  geometry.) 
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ABSTRACT 

This  is  si  discussion  of  progress  on  some  continuing 
efforts  to  understand  and  predict  the  influence  of  shock 
waves  on  turbulent  fields.  Observations  show  that  the 
average  turbulent  energy  and  Reynolds  stress  increase 
following  shock  interaction.  Strong  support  tor  ths  peak 
stress  empMcatton  is  provided  by  the  sharply  accelerated 
random  strain  field  appearing  immediately  downstream  of 
the  shock.  This  is  succeeded  by  a  gradual  decay  from 
peak  amplification  to  a  new,  more  modestly  amplified, 
near-equitibrium  level  that  persists  a  substantial  distance 
behind  the  shook.  Practical  interest  is  focused  on  the  the 
alterations  to  and  influences  of  post-shock  turbulent 
kinetic  energy,  transport  component  mixing,  wal  shear, 
and  heat  transfer.  This  work  addresses  some 
fundamental  questions  that  remain  about  the  dynamics  of 
the  amplification  process.  Initially,  viscoelastic 
turbulence  response  parameters  are  developed  from 
Monte-Cario  shock  structure  computations  to  modal  the 
localized  shock  front  interaction.  In  the  post-shock  decay 
region,  compressible  two  dimensional  large-eddy 
simulations  (LES)  are  applied. 


INTRODUCTION 

When  e  shock  moves  in  and  interacts  with  an  upstream 
turbulent  field,  the  intensity  of  ths  fluctuating  components 
(velocity,  density,  temperature)  increases  after  the 
encounter.  Turbulent  component  mixing,  transport,  and 
diffusion  are  correspondingly  enhanced. 

Basic  understanding  of  the  shock-turbulence 
interaction  dynamics  is  required  for  acuratety  estimating 
the  influence  of  the  shock-enhanced  turbulsnce  on  the 
design  of  components  exposed  to  supersonic  turbulent 
flow.  It  is  required  tor  the  development  of  useful  design- 
support  predictive  models  and  for  the  effective  design  of 
confirming  experiments. 

Considerations  about  the  physics  of  shock  turbulence 
interactions  have  a  high  priority  in  supersonic/hypersonic 
aerodynamic  flow  field  design  analysis,  exposed  surface 
materials  selection,  fabrication,  and  performance 
analysis,  as  well  as  supersonic  combustion  process  and 
design  analysis.  Shock  turbulence  interaction 


considerations  also  exert  considerable  influence  on  the. 
perhaps,  less  familiar  work  in  simulating  and  modeling  the 
evolution  dynamics  of  thermal,  electromagnetic,  and 
radiation  fields  encountered  in  astrophysicai.  magnetically 
confined  plasmadynamic.  and  inertial  confinement  fusion 
processes. 

Many  experimental  observations  of  the  enhancament  of 
shocked  turbulent  field  intensity  during  shock  turbulence 
interactions  ere  reported  tor  jets,  wakes,  free  shear  layers 
and  boundary  layer  flows. 

Some  of  the  variety  of  experimental  configurations  used 
in  association  with  the  present  work  include 
measurements  from  shock  tube  experiments  in  the 
reflected  shock  phase  (Trofiar  and  Duffy.  1S85;  Kefier  and 
Merzkirch,  1990),  shock  tube,  incident  shock  phase 
measurements  of  the  interaction  with  a  pre-existing 
randomly  varying  density  field  (Hessefink  and  Stunevant, 
1988)  and  supersonic  compression  corner  shock 
boundary  layer  interaction  experiments  (Smits  and  Muck, 
1987). 

While  qualitative  agreement  exists  on  the  enhancement 
of  the  turbulence  by  shock  transition,  the  level  of  peak 
enhancement,  the  persistence  of  the  enhancement 
downstream,  the  post  shock  energy  partition  and 
distribution,  as  weii  as  changes  to  turbulent  spectral 
structure,  length  and  time  scales  are  inconsistent, 
unresolved,  or  in  disagreement 

This  is  not  a  criticism  of  the  experiments  but  rather 
attests  to  the  basic  difficulty  inherent  in  extracting 
turbulence  information  from  compressible  supersonic 
flows.  This  type  of  experiment  imposes  severe 
requirements  on  diagnostic  technique  development  and 
the  ingenuity  of  the  experimentalist  in  designing  the 
experiments  and  interpreting  the  results.  One  is  faced 
with  the  formidable  task  of  rapidly  sampling  a  sufficient 
number  of  highly  resolved  statistical  realizations  from  a 
supersonic  background  flow  in  order  to  be  ab.«  to  interpret 
and  deduce  useful  information  from  as  many  as  three 
uncorrelated  statistical  fields. 

The  collaboration  between  those  involved  in  physical 
experiments  with  those  involved  in  numerical  simulations 
is  becoming  more  common  in  practice.  The  present  work 
on  targe-eddy  simulations  responds  to  this  theme. 
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The  computations  are  developed  in  tha  spirit  of  Bird's  - 
(1967)  one-dimensional  shock  structurs/transport 
simulations,  but  with  an  additional  (transvarsa)  dimansion 
addad  for  angular  momentum  exchange,  curvature  and 
vortidty  development. 

The  statistical  'reservoir  states*  for  the  ensemble  of 
interactive  encounter  values  describing  upstream  and 
downstream  conditions  had  to  be  increased  appropriate  to 
the  additional  degrees  of  freedom.  Pre-shocked  and  post- 
shocked  states  and  the  evolution  of  the  distribution  of 
states  over  the  shockwave  interaction  region  are 
simulated  in  a  discretized  macroscopic  analogy  to  the 
asymptotic  stats  conditions  derived  in  Mod-Smith's  (1951) 
theoretical  expansion  procedure  for  kinetic  theory  shock 
structure  calculations. 

Additional  motivation  for  these  simulations  was  provided 
by  the  theoretical  observation  of  the  viscoelastic 
response  evident  for  an  isotropic  turbulence  field  on 
imposition  of  a  sudden  variation  in  mean  strain.  This 
finding  was  suggested  in  analyses  of  the  evolution  of 
turbulent  field  structure  using  direct  interaction 
approximation  (DIA)  computations  (Kraichnan,  1967). 

Figure  3.  shows  the  filtered,  ensemble-averaged 
probability  density  distributions  of  fluctuating  velocity 
(u*{),  internal  energy  (E*).  and  mass  density  (p*) 
fluctuations  computed  through  a  distance  swept  by  an 
unsteadily  deforming  shock  front  during  interaction  with 
turbulence.  Here  we  are  concerned  with  ensemble 
averaged  peaks  (designated  with  a  double  asterisk)  of  the 
distribution  of  the  turbulence  pressure,  p**»  <p**xq2>, 
density  fluctuations,  <p**>,  and,  by  impication.  velocity 
fluctuations  since  <q2>  -  1/2<uj"2  +  Uj”2>.  The 
brackets,  o.  designate  ensemble  averages  taken  over  5 
averaged  realization  distribution  profies. 

The  shock  deforms  and  vibrates  about  its  mean  position 
(in  a  coordinate  system  moving  with  a  Mach  3  shock  front) 
as  shown  in  the  shaded  silhouette  trace  on  the  right  side 
of  the  figure.  This  trace  is  the  spatial  displacement  swept 
by  the  shock  from  fluctuations  over  the  total  Monte  Carlo 
simulation  period.  Details  of  the  procedure  and  the  results 
are  reported  in  Buckingham  (1989, 1990). 

In  preliminary  tests,  the  Monte-Cario  simulations  were 
applied  for  comparison  to  the  Shigemi,  Koyama,  and 
Alhara  (1976)  unsteady  shock  resonance  tube 
experiments.  Results  of  simulations  on  the  shock 
oscillation  frequency  as  well  as  the  entropy  defect  in 
comparison  to  stationary  Hugoniot  conditions  (or 
approach  to  adiabatidty  by  unsteady  shock  waves)  were 
compared  with  the  resonance  tube  experiments.  The 
Monte  Carlo  results,  which  indude  both  the  influence  of 
random  deformation  of  the  from  as  well  as  random  motion 
about  the  mean  shock  position  somewhat  overpredict  the 
energy  transferred  to  turbulence,  but  the  predicted  values 
and  qualitative  behavior  (the  apparent  rate  of  decay  of  the 
rate  of  transfer  of  directed  shock  energy  to  turbulence 
energy  with  increasing  Mach  No.)  appears  to  be 
reasonable  representations  of  the  experimental 
results. This,  however,  is  a  candidate  for  more  systematic 
future  investigations  using  a  broader  experimental  data 
base,  i  available. 


Table  I  summarizes  some  of  the  dimensional  scale 
results  of  the  Monte-Carlo  simulations.  Here  the 
diminishing  unsteady  shock  interaction  motion  range  is 
traced  from  Mach  2  to  Mach  6  in  the  Monte-Carlo 
simulations.  This  is  designated  D  in  the  table.  The  values 
are  ratioed  to  a  unit  unsteady  displacement  trace  at  Mach 
2.  The  correlation  length  scales  .A,  of  the  principle  eddy 
distortions  of  the  shock  from  are  given  in  the  second 
column,  again  ratioed  to  the  unit  value  assigned  to  the 
trace  of  the  overall  shock  motion  at  Mach  2.  This  length 
scale,  A.  is  an  important  parameter  in  the  viscoelastic 
response  model  frequency  discussed  later.  The  ensemble 
averaged  amplification  of  the  preshocked  turbulence 
streamwise  kinetic  energy  component  (initiated  at  4% 
relative  intensity)  is  given  in  the  third  column. 

We  move  now  to  the  general  solution  of  the  resolved  grid 
scale  motions  and  modifications  used  for  the  Smagorinsky 
subgrid  scale  eddy  viscosity  model.  The  solution  is 
approximated  on  the  discretized,  time  explicit,  two- 
dimensional  equations  of  compressible  viscous  flow, 
where  averaging  is  mass-weighted  (Favre,  1983)  over  the 
Smagorinsky  modeled  subgrid  scale  motions  represented 
by  a  total,  shear,  buoyancy  (bulk)  and  dilatation 
contribution  to  the  eddy  viscosity,  vt. 

The  equations  are  written  in  conservation  form  prior  to 
the  discretization  defining  spatial  mesh  distribution  and 
explicit  time  evolution.  For  simplicity,  this  description  will 
be  confined  to  a  regular  Cartesian  coordinate  space  in 
contrast  to  the  variety  of  coordinate  transformations  used 
in  applications  involving  complex  geometries.  The 
variable  labels  are:  coordinate  directions.  x$,  velocity 
components  and  their  derivatives  with  respect  to  the 
coordinate  directions,  U|,  Ujj.  mass  density,  p,  and 
specific  total  energy,  E  -  e  +  uk^  /  2.  Here  e  is  the  specific 
internal  energy.  A  calorically  corrected  polytropic 
equation  of  state  of  the  gas  used  for  these  trials;  with 
static  pressure,  p»  (y-l)pe,  temperature,  8  •  e/Cy,  where 
cv  is  corrected  for  bound  stale  excitation.  The  averaged 
mass,  momentum,  and  energy  conservation  equations  for 
the  resolved  scales  are: 

P.t  ♦  (puk).k  -  0.  (1) 

(puj),t+  <pujuk).k  -  -p.j+  (pvt  Sjk).k  «■  ptejj^j)*.  (2) 
(pE).t  +  (pEuk),k  .  ((pv(/n*)E^).k  +  pvtSikUj,k 

*  P°k*  *  P<£ijk$j)*kui-  (3) 

The  tensor  summation  convention  holds  in  the  foregoing 
equations,  is  the  standard  alternating  tensor  and  Sjk  is 
the  Kronecker  delta  tensor. 

The  resolved  scale  deviatoric  strain  rate  tensor  on 
which  the  local  shear  production  depends  is  given  by, 

S*  -  Uj*  +  u^j  -  2/3uj.j8jj 
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Away  from  the  shock,  at  modest  finite  distance,  E,  the 
form  of  the  eddy  frequency  given  by  equation  (7)  obtains 
directly. 

The  generalized  eddy  viscosity  is  given  by, 

v{  -  (CsX)2f.  (10) 

Results 

Predicted  LES  from  Buckingham  (1991,  symbol  B)  and 
from  Rotman  (1991.  symbol  R)  shock  amplification  of  the 
turbulent  kinetic  energy  over  the  total  grid  is  illustrated  in 
Figures  4  and  5  in  comparison  to  experimental  values  from 
Keller  and  Merzkirch  (1990,  symbol  K  &  M).  Figure  4 
shows  pre-ehocked  values  of  power  spectral  density  while 
Figure  S  shows  the  post-shocked  values.  The  power 
spectral  densities  which  are  different  for  the  experiments 
and  the  simulations  are  rescaled  and  renormalized 
individually  to  the  highest  decade  obtained  in  their  post¬ 
shock  value.  Hence,  spectral  shape  and  relative  change 
may  be  obtained  from  these  figures  instead  of  absolute 
quantitative  values. 

Here  the  Rotman  results  show  some  evidence  of  long 
period  computational  dissipation  in  the  absence  of  a 
subgrid  scale  model  while  predicting  an  increase  in  the 
overall  grid  spectral  energy  content  of  about  20%  following 
shock  transition.  The  Buckingham  results  predict  about 
an  85%  increase  in  the  kinetic  energy  while  the 
experimental  values  at  these  flow  conditions  reflect  a 
measured  increase  of  about  150%  where  limited 
measurement  site  distributions  may  be  reflected  in  this 
disparity.  However,  these  results  seem  to  confirm  the 
utlity  of  a  formal  procedure  for  modeling  the  influence  of 
the  non-resolvable  subgrid  scales  even  when  using 
advanced,  high  resolution  multi-grid  shock  capturing 
procedures  for  the  resolved  scales  in  LES. 

The  3  compression  comer  cases  illustrated  in  Figures  6, 
7,  and  8  serve  to  illustrate  the  need  for  recent  efforts  on 
implementing  a  more  general  subgrid  scale  model  for  near 
wall  anisotropic  compressible  behavior. 

In  Figure  6  the  effects  of  boundary  layer  diffraction  and 
near  wall  dissipation  on  the  sharply  inclined  oblique 
shockwave,  supported  by  an  8  degree  half  angle 
cc  T’Dression  comer,  are  apparent  The  actual  outermost 
st  «amiine  kinetic  energy  profile  (at  y*-0.6)  taken  from 
the  experiments  of  Smtts  and  Muck  (1987)  shows  a  much 
more  modest  amplification  than  the  the  results  predicted 
by  the  LES  scheme.  Furthermore,  the  near-  wall 
influences  at  y*<  0.6  could  not  be  successfully  modeled. 
The  outermost  streamline  is  the  only  one  in  which  the 
results  are  comparable  in  the  neighborhood  of  the  waJL 
Results  for  a  16  degree  half  angle  compression  comer, 
are  shown  in  Figure  7.  Somewhat  more  satisfactory 
agreement  is  obtained  (in  part,  fortuitously)  because  of 
the  emergence  to  a  less  shallow  inclination  of  the  shock  in 
the  shock  boundary  layer  interactive  region.  Again, 
results  for  the  lower  lying  streamlines  indicate  the  need  of 
improvement  in  the  dissipative  near-wall  region. 

Figure  8.  indicates  that  for  even  steeper  shocks  the 
outermost  streamline  results  are  satisfactorily  reasonable 
in  comparison  to  experiment  It  should  also  be  noted  that 
these  current  results  represent  ensemble  averages  over 


at  least  five  computational  runs  as  opposed  to  the  very 
inadequate  results  of  single  realizations  previously 
shown.  The  confirms  a  suggestion  made  by  Leith  (1991) 
on  analyzing  the  results  of  a  properly  set  LES  resolved 
scale  data  set  as  stochastic,  with  interpretations  based 
on  ensemble  averages  of  a  repeated  set  of  realizations. 

SUMMARY  AND  PROJECTIONS 

More  generally  effective  compressible  subgrid  scale 
models  may  be  obtained  by  including  the  intensity  of  the 
density  fluctuations  as  a  parameter,  together  with  a 
turbulent  Mach  number  and  including  the  influence  of  the 
compressible  dilatation  explicitly  as  pointed  out  in  the 
analysis  of  the  subgrid  scale  closure  process  by  Speziale. 
Eriebacher,  Zang  and  Hussaini  (1988).  In  this  work 
advantage  is  taken  of  the  appropriate  use  of  low  Reynolds 
Number  direct  numerical  simulations  to  illustrate  some 
crucial  features  of  the  wave  dominated  turbulent  transport 
region  in  comparison  to  the  rotationaliy  dominated  region. 
Implementation  and  later  encouraging  results  in  the 
evolution  of  the  LES  schema  is  provided  in  Eriebacher, 
Speziale.  Hussaini  and  Zang  (1990). 

One  of  the  several,  more  universal  subgrid  model 
developments  is  a  two  band  grid  to  subgrid  scale  overlap 
model  where  an  algebraic  identity  provides  a  procedure  for 
making  use  of  the  resolved  scale  spectrum  in  dynamically 
altering  subgrid  scale  eddy  viscosity  coefficients 
(Germono,  Pnmelli,  Mom.  and  Cabot  1991).  This  has 
been  extended  to  compressible  flows  more  recently  by 
Mom.  Squires,  Cabot  and  Lee  (1991).  Implementation 
and  test  are  delayed  by  properly  adjusting  average 
frequencies  along  at  least  one  honrogeneous  dimensional 
direction.  In  the  present  compression  comer  experiments 
this  is  represented  by  the  spanwise  direction.  This  work  is 
in  progress,  but  is  inconclusive  at  this  stage. 

It  is  suggested  that  the  shockwave  turbulence 
interaction  problem  is  a  prototypical  physical  process  that 
is  imperfectly  understood  and  yet  dominates  many  typical 
supersonic  and  hypersonic  flow  investigations.  It  is  also  a 
daunting  numerical  problem  combining  the  special 
problems  of  resolving  shockwaves  and  turbulence.  This 
discussion  emphasizes  current  progress  in  developing  all 
Reynolds  no.  compressible  LES  procedures  to  assist  us  in 
our  investigation  of  the  underlying  physical  mechanisms 
that  govern  the  interaction.  The  discussion  here  focuses 
on  the  progress  to  date  in  tne  search  for  effective  subgrid 
scale  closure  procedures  and  numerical  algorithm 
development.  The  search  is  dearly  far  from  over. 
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FIGURE  1.  ROTMAN  (1991)  ADAPTIVE  MESH  GODUNOV 
LES  DENSITY  GRADIENT  ILLUSTRATION  OF  M  -  1.5 
SHOCK  PASSAGE  THROUGH  4.7%  RELATIVE 
INTENSITY  DENSITY  FLUCTUATIONS. 


FIGURE  3.  FILTERED  AVERAGE  PROBABILITY  DENSITY 
DISTRIBUTIONS  OF  FLUCTUATIONS  IN  INTERNAL 
ENERGY  (E).  DENSITY  (p)  AND  VELOCITY  (U). 


FIGURE  2.  DARK  FELD  SHADOWGRAPH  OF  LASER¬ 
GENERATED  MACH  100  SHOCK  OVER  RUNNING  AND 
DISTORTED  BY  PASSAGE  THROUGH  TURBULENCE. 


FIGURE  4.  NORMALIZED  POWER  SPECTRAL  DENSITY 
OF  UPSTREAM  TURBULENCE  COMPARING  COMPUTED 
RESULTS  FROM  :  ROTMAN  (1991),  E1(R);  CURRENT 
RESULTS.  E2(B);  AND  EXPERIMENTS  OF  KELLER  AND 
MERZKIRCH  (1990),  E3(K&M). 
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ABSTRACT 

Presented  are  recent  computational  results  in  continuing  work  on  the 
description  and  prediction  of  the  influence  of  shock  wave  interactions  on 
turbulence.  Depending  on  the  strength  of  the  shock  and  the  intensity  of  the 
pre-shocked  turbulence,  interaction  may  significantly  increase  turbulent 
energy  and,  in  some  situations,  severely  distort  the  shockwave.  Behind  the 
shock  front,  the  amplification  decays  to  a  new,  mote  modestly  amplified  state, 
relative  to  the  pre-shocked  level.  Practical  consequences  include  the  influence 
of  this  shock  turbulence  enhancement  on  post-shock  eddy  transport, 
component  mixing,  diffusion,  surface  shear  and  heat  transfer.  In  the  present 
work,  we  combine  a  viscoelastic  response  model,  developed  from 
independent  Monte-Carlo  shock  interaction  simulations  with  compressible, 
large-eddy  simulations  (LES).  Results  are  compared  and  interpretations 
developed  with  the  aid  of  experiments,  including  some  remarkable  new  data 
for  very  strong,  hypervelocity  shocks  interacting  with  intense  turbulence. 


1.  INTRODUCTION 

When  a  shock  interacts  with  a  pre-existing  turbulent  field,  a  substantial 
enhancement  of  the  turbulent  intensity  and  a  simultaneous  distortion  and 
realignment  of  the  turbulent  strain  rate  field  and  resulting  turbulent  stress  field 
occurs.  Component  mixing,  and  transport  are  often  significantly  enhanced. 

Combinations  of  low  upstream  turbulent  intensity  and  strong  shocks 
bring  about  the  most  pronounced  changes  in  the  level  of  turbulent  intensity. 
Contrastingly,  for  weak  shocks  and  intense  upstream  turbulence,  shock 
distortions,  and  localized  shock  surface  modons  are  most  evident.  However, 
very  recent  experimental  evidence  suggests  that  for  intense  upstream 
turbulence,  significant  shock  distonions  may  develop  in  very  strong, 
hypervelocity  shock  fields  [1]. 
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Accurate  description  of  the  shock-turbulence  interactions  is  vital  for 
development  of  useful  predictions  and  predictive  models  of  their  influence  on: 
supersonic/hypersonic  flow  field  analysis,  design,  and  materials  selection, 
interior  supersonic  exterior  flow  field  analysis,  design,  materials  selection, 
interior  supersonic  combustion  analysis  and  combustion  chamber  design.  It  is 
also  crucial  for  accurately  predicting  the  development  and  evolution  of  flow 
Held  generated  thermal  and  electromagnetic  radiation  fields  which  are 
important  considerations  in  supersonic/hypersonic  configuration  design 
analysis. 

The  shock  tube  experiments  of  Hesselink  and  Sturtevant  [2]  illustrate 
that  substantial  shock  front  deformation  may  occur  when  a  weak  shockwave 
traverses  a  relatively  intense  initial  turbulent  field.  The  interactions  may  even 
eventually  lead  to  the  apparent  shock  front  break  up  of  weak  shocks  into 
separate  waves.  These  experiments  are  the  focus  of  Rotman's  [3]  LES 
investigations  discussed  later.  Emphasis  was  placed  on  capturing  and 
resolving,  as  well  as  possible,  the  distortions  of  the  shock  front  during 
interactive  passage  through  turbulence.  The  multi-grid  Euler  procedure  of 
Colella  and  Glaz  [4]  was  adapted  for  these  high-resolution  shock  interaction 
results  using  an  initiated  random  velocity  and  density  field. 

In  contrast,  it  is  a  commonly  held  view  that  if  the  interactive  shock  Mach 
number  is  in  the  hypersonic  range  (shock  Mach  numbers  equal  to  or  greater 
than  about  6),  rapid  reinforcement  of  shock  strength  at  the  front  by  non-linear 
pressure  Held  interactions  and  consequent  self-healing  of  shock  front 
geometric  irregularity  act  to  suppress  front  distortion  and  also  suppress  the 
influence  of  the  shock  interactions  on  turbulence.  In  the  present  work,  recent 
experiments,  numerical  large  eddy  simulations  (LES),  and  analysis  are 
presented  which  contradict  this  notion.  Pronounced  shock  front  distortion  and 
substantial  interactive  modification  of  the  turbulence  seems  to  persist  even  for 
shocks  propagating  at  shock  Mach  numbers  well  in  excess  cf  100. 

A  primary  goal  of  the  present  work  is  to  provide  a  supplementary,  high 
precision,  numerical  simulation  tool  for  augmenting  the  sparse  experimental 
evidence  on  shock  turbulence  interactions  at  these  elevated  Mach  Numbers. 
Some  commonly  held  theoretical  views  (such  as  those  described  in  the 
previous  paragraph)  appear  to  be,  at  least,  inadequately  descriptive  of  the 
limited  experimental  data  and,  more  significantly,  may  be  "potentially 
misleading.  Results  of  systematic  numerical  simulations  for  selected 
geometries  and  flow  condidons  may  substantially  improve  this  situation. 


2.  FORMULATION  AND  NUMERICAL  PROCEDURE 

The  2D  subgrid-scale  model  used  for  most  of  the  present  work  is  the 
simple  algebraic  Smagorinsky  model  [5]  with  assumed  isotropic 
production/dissipation  balance,  but  adjusted  at  the  sub-grid  level  for 
compressibility  sources  and  dissipation.  Point-to-point  seeding  of  directly 
simulated  space  and  time  white  random  perturbational  acceleration  provides 
the  initial  disturbance.  Evolution  to  the  desired  average  initial  turbulence  field 
is  computed  over  a  predetermined  initial  trial  computational  phase.  This 
process  is  repeated,  as  necessary,  with  slight  alterations  to  the  ensemble 
statistics  until  the  prescribed  initial  preshocked  turbulent  state  develops. 


28 


Favre  [9]  mass  averaging  is  included  in  the  discretized  grid  scale  Navier 
Stokes  equations  which  are  written  in  conservation  form.  The  general  LES 
formulation  follows  directly  from  the  model  developments  outlined  for 
studies  of  the  2D  shear  layer  in  Leith  [6]  and  for  the  compressible  3D 
turbulence  massively  parallel  processor  simulations  of  Leith  [7]. 

A  dynamic  viscoelastic  response  model  for  the  turbulent  field  changes 
on  shock  interaction  couples  directly  to  the  shock-free  Smagohnsky  subgrid 
model  through  buoyancy  source,  divergence  work  and  dilatational  dissipation 
relations.  Coupling  occurs  at  the  sub-grid  scale  characteristic  frequency  level. 
Use  is  made  of  relaxation  and  instantaneous  modulus  parameters  developed  in 
a  previous  phase  of  the  present  work  through  application  of  the  results  of 
direct  Monte-Carlo  shock  structure  simulations  [8,9]. 

Perturbational  variables  developed  from  the  Monte-Carlo  simulations 
include  ensemble-averaged  probability  density  distributions  of  fluctuating 
velocity  (u*j),  internal  energy  (E*),  and  mass  density  (p*).  Die  distributions 
are  computed  over  a  spatial  range  of  statistical  shock  turbulence  interaction 
equivalent  to  the  physical  distance  swept  by  the  deforming  shock  front  during 
interaction  with  turbulence.  For  LES  model  application,  use  is  made  of 
specific  ensemble  averaged  peak  values,  as  well  as  correlations,  and  integral 
scales.  The  peak  values  are  designated  with  a  double  asterisk.  Consideration 
is  given  to  peak  turbulent  pressure,  p**=  <p**xq2>;  peak  fluctuations  in 
mass  density,  <p**>;  and,  by  implication,  peak  component  velocity 
fluctuations,  where  <q2>  =  l/2<uj**2  +  uj**2>.  Here  the  brackets,  <>, 
designate  ensemble  averages  usually  formed  over  5  realization  distributions. 

The  Monte-Carlo  simulations  also  yield  required  correlation  length 
scales  of  the  principle  shock  front  eddy  distortions,  A.  This  length  scale  is  an 
important  parameter  in  the  viscoelastic  response  model  frequency  discussed 
later. 


We  move  now  to  description  of  the  general  procedure.  The  LES  grid 
scale  (resolved)  motions  are  developed  from  discrete  approximations  to  the 
time  dependent,  two-dimensional  equations  of  compressible  viscous  flow. 
Favre  [9]  mass-weighted  averaging  applies  to  the  resolved  motion  scales 
down  to  a  dimension  of  the  order  of  the  mesh  spacing.  The  influence  of  the 
unresolved  scales  are  modeled  as  single  turbulence  wave-length  band  grid 
scale  averages  or  new  two  band  scale  overlap  procedures  at  the  scale  division 
separating  computed  from  unresolved  scales  of  motion. 

The  latter  yields  a  procedure  for  dynamic  modifications  to  the 
compressible  Smagorinsky  sub-grid  scale  model  that  adjusts  modeled  near 
field  behavior  for  wall  influences  and  provides  additional  degrees  of  freedom 
for  channeling  the  influence  of  the  grid  scale  to  the  modeled  subgrid  scale 
motions.  Test  applications  of  the  two  band  model  are  currently  in  progress. 
[10].  Here  the  subgrid  scale  compressibility  is  modeled  by  shear,  buoyancy 
(bulk),  and  dilatation  contributions  to  the  eddy  viscosity,  vt. 

For  simplicity,  this  description  will  be  confined  to  a  regular  Cartesian 
coordinate  space,  in  contrast  to  the  variety  of  coordinate  transformations  used 
in  applications  involving  complex  geometries.  The  variable  notation  used 
includes:  coordinate  directions,  xj;  velocity  components  and  their  derivatives 
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with  respect  to  the  coordinate  directions,  uj,  uj.  j;  mass  density,  p;  and 
specific  total  energy,  E  =  e  +  u^2  /  2.  Here  e  is  the  specific  internal  energy. 

A  calorically  corrected  quasi-ideal  polytropic  equadon  of  state  of  the  gas 
is  used  for  the  cunent  tests,  in  which  the  stadc  pressure  is  given  simply  by, 
P  =  (Y-l)pe.  The  caloric  correction  applies  to  the  temperature,  0  =  e/cv, 
where  cv  is  corrected  for  bound  state  or  continuum  excitation. 

The  averaged  mass,  momentum,  and  energy  conservation  equations  for 
the  resolved  scales  are  written: 

P-t  +  (puk).k  =  0,  (1) 


(puj).t  +  (pUjUk).k  =  -  P.i  +  (pvt  Sjk),k  +  p(£ijk<J>j).k.  (2) 

(pE),t  +  (pEuk),k  *  [(pvt/n*)E.k3.k  +  pvtSfcUi.k 


-  puk,k  •  p(ejjk<t>j).kUi-  (3) 

The  tensor  summation  convention  holds  in  the  foregoing  equations,  £jjk 
is  the  standard  alternating  tensor.  Also,  Sy  denotes  the  Kronecker  delta 
tensor. 

The  resolved  scale  deviatoric  strain  rate  tensor  on  which  the  local  shear 
production  depends  is  given  by, 

Sik  *  Uj.k  +  uk,j  •  2/3uj,j5jj.  (4) 

In  some  of  the  present  computations  tests  are  underway  on  stochastic 
backscatter  influences.  [101  These  influence,  to  lowest  order,  the  smallest 
resolved  scales  (marginally  larger  than  the  resolution  scale  of  the  grid). 
Definitive  comments  await  further  analysis. 

However,  for  completeness,  we  introduce  the  form  of  the  stochastic 
modeling  used  for  testing  these  influences  at  this  time.  The  influences  are 
explicitly  modeled  using  the  space  and  time  white  random  tensor  acceleration 
potential  developed  by  Leith  [6,7], 

*k  =  Cb(*8t)3/2(X/5t)2gK,  (5) 

This  is  applied  at  each  grid  point  and  at  each  explicit  time  step  in  the 
calculation,  where  the  characteristic  resolution  length  scale,  X  ,  is  taken  to  be 
twice  the  grid  scale.  Ax.  This  acts  as  a  smoothing  filter  to  partially 
compensate  for  coarse  grain  random  errors  that  develop  in  propagating 
random  disturbances  at  exactly  the  resolution  scale.  The  explicit  time  step  of 
the  calculation  is  represented  by  5t  The  vector  term,  gK,  components  are  unit 
gaussian  random  numbers,  each  drawn  from  a  population  with  zero  mean  and 
unit  variance.  The  generalized  eddy  frequency,  f,  and  the  constant,  Cu  are 
identified  in  the  subsequent  Smagorinsky  subgrid  scale  model  discussion. 


30 


An  implicit  assumption  is  that  the  velocity  field  in  a  test  problem 
contains  an  indefinitely  large  number  of  scales  of  motion  (a  unit  Reynolds 
number  of  the  order  of  perhaps  105  or  greater),  so  that  the  classical 
Kolmogorov  inertial  range  separation  between  production  and  dissipation 
emerges  with  a  universal  energy  spectrum,  E(k)  =  aE^tr5/3.  This  provides 
the  necessary  basis  for  the  dimensional  analysis  scaling  and  arguments  used 
for  estimating  the  energy  dissipation,  E.  In  the  following  outline  K  represents 
the  total  specific  turbulent  energy,  E(k),  integrated  over  the  entirety  of  wave 
number,  K,  space. 

In  the  simple  eddy  viscosity  concept,  the  local  eddy  stress  tensor  is 
directly  related  to  the  averaged  rate  of  strain  tensor  field.  In  large  eddy 
simuladons,  this  averaged  strain  rate  field  consists  of  the  explicitly  resolved 
scales  of  motion. 


Tjj  =  vtsij-  (6) 

and  the  shear  production  of  the  turbulent  kinetic  energy  is  given  by  the  work 
in  producing  the  eddy  viscous  stress, 

K,t  -  TjjUj.j  =  vtS2,  (7) 

where  the  mean  resolved  strain  rate,  S,  is  defined  so  that, 

S  -  (SijU|.j)1/2  =  <eijsi/2>1/2  (8) 

For  compressible  flow,  the  incompressible  Smagorinsky  balance 
condidon  between  production  and  dissipation  of  the  energy  must  be  modified 
for  the  compressible  buoyancy  source,  B,  of  kinetic  energy  (which  reaches  a 
peak  at  the  passage  of  a  shock  wave  ),  as  well  as  the  dilatational  turbulence 
energy  lost  or  dissipated  by  turbulent  pressure  work  on  local  velocity 
divergence,  KD.  Here  D  *  uj.j,  while  the  buoyancy  source  term,  B,  is 
computed  from  the  product  of  buoyancy  times  acceleration, 

B  s  -  Sc'1  (p,j/p)(p,j/p).  (9) 

Here  Sc  is  a  Schmidt  number  which  was  assigned  a  value  of  0.7  as  in 
Leith  [6,7].  B  is  positive  or  zero.  In  shock  interaction  simulations,  B 
represents  the  computed  transient  positive  definite  viscoelastic  modulus  based 
on  instantaneous  state  and  relaxation  parameters  developed  from  direct 
Monte-Carlo  shock  interaction  simulations  [8,9].  The  compressibility 
modified  Smagorinsky  balance  between  production  and  dissipation  is  written, 

vt(S2  +  B)  -  2/3KD  =  e,  (10) 

from  which  the  generalized  eddy  frequency  away  from  the  shock  may  be 
evaluated, 


f  «  (S2+  B  +  Cd2D2)1/2  -  CdD.  (11) 
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As  observed  by  Leith  [7],  this  is  equivalent  to  the  familiar  Von- 
Neumann  Richtmyer  shock  damping  dissipation  term  when  the  dimensionless 
coefficient,  Cd,  is  set  equal  to  10. 

At  the  shock,  we  define  the  characteristic  viscoelastic  interaction 
response  frequency  from  the  parameters  identified  earlier, 

G  =  [  B  + 1  g  -  B  I  exp  -  (E/A)n],  (12) 

g  *  p*7(<p**>2A2). 

In  the  exponential  relaxation  term,  £,  represents  the  stream  wise  distance 
behind  the  mean  shock  front  position  and  the  exponent,  n,  takes  the  value  2.5 
in  our  present  studies,  based  on  the  results  of  the  stochasdc  shock  interaction 
calculations  discussed  earlier.  Near  the  shock,  the  generalized  eddy  frequency 
becomes. 


1  =  (S2  +  G  +  Cd2D2)1/2  .  cdD  (13  ) 

Away  from  the  shock,  at  modest  finite  distance,  £,  the  form  of  the  eddy 
frequency  given  by  equation  (7)  obtains  directly. 

The  generalized  eddy  viscosity  is  given  by, 

vt=(CsX)2f.  (14) 


3.  RESULTS  AND  COMPARISONS  WITH  EXPERIMENTS 

We  initiate  this  pan  of  the  discussion  with  results  of  simulations  of 
shock  boundary  layer  interactions  on  three  two-dimensional  wedge 
compression  comers  at  a  nominal  Mach  number  of  2.7,  in  air.  The  three 
compression  comer  half  angles  are  8,  16  and  20  degrees.  Detailed 
experimental  average  turbulence  profiles,  wall  distributions  and  statistical 
structure  boundary  layer  information  are  taken  from  the  experiments  of  Smits 
and  Muck  [12].  In  the  interest  of  conserving  space  in  this  paper,  one  example, 
that  of  the  16  degree  half  angle  test  case,  will  be  used  for  our  discussions.  In 
the  experiments  and  the  simulations  an  oblique  shock  forms  at  the  upstream 
compression  comer.  The  foot  of  the  shock  is  immersed  in  a  turbulent  wall 
boundary  layer  layer,  with  which  it  interacts.  Measured  values  of  the 
turbulent  kinetic  energy  distribution  downstream  of  shock  interaction  ratioed 
to  the  measured  upstream  values  are  used  as  a  basis  of  comparison  and  test  of 
the  LES  plus  shock  interaction  model. 

In  the  experiments  and  the  LES  trials  the  distributions  of  turbulence 
kinetic  energy  are  traced  along  three  streamline  paths  displaced  at  y*  =  0.2, 
0.4  and  0.6  from  the  wall  surface,  where  y*  is  the  ratio  of  the  normal  distance 
from  the  wall  to  the  total  boundary  layer  thickness  at  the  depicted  fluid 
element  transit  time  following  shock  interaction.  Figure  1.  shows  the 
experimental  vs.  LES  traces.  It  indicates  that  the  computed  turbulence 
distributions  are  reasonably  comparable  for  the  outermost  streamline.  This 
particular  displacement  is  the  most  removed  from  the  near  wall  influences. 
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Shock  -  boundary  layer  TKE  amplification 
g  16  degree  wedge  compression  corner 


Time  after  shock  passage,  microsecs. 


Figure  1.  Comparison  of  current  LES  results  with  experiment  [ 12 }. 
Amplification  of  turbulent  kinetic  energy  by  shock  boundary  layer 
interaction  in  a  16  degree  compression  corner  at  Mach  2.7. 


The  result  emphasizes  the  effectiveness  of  the  viscoelastic  shock  response 
model  in  tracing  the  amplification  and  relaxation  of  the  shocked  turbulent 
field.  At  the  same  time,  the  results  indicate  the  ineffectiveness  of  the  implicit 
isotropic,  statistically  homogeneous,  subgrid  scale  model  used  here  when  in 
the  neighborhood  of  the  strongly  inhomogeneous,  anisotropic  wall  influences. 

Our  focus,  it  must  be  remembered,  is  on  developing  an  effective 
technique  which  automatically  adjusts  to  the  influence  of  both  shock  wave 
interactions  and  complicated  boundary  or  interface  geometries  on  turbulence. 
To  this  end.  we  implement  and  test  a  dynamically  adjusted  sub- grid  scale 
model  introduced  by  Germano,  et  al  [13]  and  extended  for  compressibility  by 
Moin  cl  al.  [14].  These  later  developments  are  currently  being  evaluated  as  a 
means  to  automatically  resolve  features  both  in  the  immediate  neighborhood 
of  the  shock  front  and  the  wall  influence  regions  with  concurrent  attention  to 
the  influence  of  the  unresolved  scale  stochastic  backscatter  on  the  resolved 
grid  scale  motions.  [10], 

Our  current  LES  predictions  [9, 10]  (symbol  B)  are  displayed  in  Figs.  2 
and  3  together  with  those  from  the  specialized,  precise  shock  resolution 
studies  of  Rotman  [3]  (symbol  R)  which  produced  satisfactory  comparisons  to 
low  Mach  number  shock  tube  experimental  interaction  results  [3]. 
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Figure  2  shows  pre-shocked  spectral  distributions  of  power  spectral 
density  while  Fig.  3  shows  the  post-shocked  spectral  distribution.  The  figures 
illustrate  that  predominant  shock  amplificadon  and  energy  transfer  occur  at 
the  low  wave  number  (explicidy  computable  grid  scale)  dynamic  range  of 
greatest  interest  in  both  our  LES  development  and  our  shock  interaction 
analyses.  For  comparison,  we  also  display  low  Mach  number  experimental 
shock  turbulence  interaction  spectral  distributions.  These  are  taken  from  the 
speckle  photographic  experimental  results  for  a  reflected  low  Mach  number 
plane  shock  passing  back  through  grid-generated  turbulence  created  in  the 
wake  of  the  incident  shock  phase  [11]  (symbol  K  &  M). 


Figure  2.  Comparison  of  pre-shocked  normalized  power  spectral 
density  distributions  from  current  LES  computations,  previous 
computations  [3]  and  shock  tube  reflected  shock  interaction 
experiments  111]  at  Mach  1 2. 


While  the  Mach  numbers,  initial  turbulent  states,  dynamic  range,  and  gas 
composition  were  somewhat  different  in  the  two  sets  of  computational 
simulations  and  the  underlying  experiments,  for  comparative  purposes  the 
results  are  rescaled  and  renormalized  to  individual  peak  values  taken 
independendy  from  each  of  the  experimental  or  simulated  trials.  Specifically, 
all  ordinate  quantities  plotted  are  power  spectral  density  values  vs  wave 
number  ratioed  to  the  peak  post-shocked  decade  value  obtained  for  each 
individual  experiment  or  simulation. 
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Figure  3.  Comparison  of  post-shocked  normalized  power  spectral 
density  distributions  from  current  LES  computations,  previous 
computations  [3]  and  shock  tube  reflected  shock  interaction 
experiment  s  [11]  at  Mach  12. 


In  Figs.  2  and  3,  the  Rotman  results  appear  to  show  some  evidence  of 
very  slight,  late  time  computational  dissipation  in  comparison  to  the  current 
results.  This  may  reflect  the  absence  of  explicit  control  by  subgrid  scale 
modeling.  The  integrated  energy  increase  through  shock  interaction  predicts 
an  increase  in  the  overall  grid  spectral  energy  content  of  about  20%,  following 
shock  transition.  This  is  in  reasonable  agreement  with  the  underlying 
experiments  simulated  [2], 

The  current  LES  results  predict  about  an  85%  increase  in  the  energy 
content  through  shock  amplification,  while  the  experimental  values  at  these 
flow  conditions  [11]  reflect  a  measured  increase  of  almost  double  this  amount. 
The  disparity  is  not  considered  crucial  because  of  the  limited  dynamic  range 
of  the  simulations  in  comparison  to  the  experiments.  Both  experimental  and 
computational  cascades  are  seen  to  be  somewhat  steeper  than  the  classical 
k~5/3  two  dimensional  spectral  decay  law,  possibly  reflecting  some  3 
dimensional  influences  in  the  experiments  and  some  uncontrolled  numerical 
dissipation  in  the  simulations.  Overall  the  behavior  appears  quite  reasonable. 
We  move  on  to  even  more  promising  results  in  the  hypervelocity  shock 
turbulence  interaction  range. 

The  most  recent  experiments  [1],  used  as  a  primary  data  base  for  our 
latest  LES  developmental  efforts,  provide  unique  and  remarkable 
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experimental  evidence  of  the  existence  of  pronounced  interactive  shock 
distortion  and  turbulent  field  alteration  even  for  a  very  strong  (Mach  number 
>  100)  spherical  shock  wave.  The  experimental  findings  appear  to  be  well 
supported  by  our  current  LES  simulations,  created  to  assist  in  the 
interpretation  and  analysis  of  the  experimental  data.  In  our  experiments,  a 
4-ns  pulse  from  the  1.06-(J.m,  1  kJ  Pharos  III  laser  at  the  Naval  Research 
Laboratory  is  focused  onto  the  surface  of  a  period-sized  piece  of  target 
material  placed  in  an  ambient  gas  of  optionally  selected  composition  and 
background  density.  The  laser  heats  the  material  to  a  few  hundred  eV, 
creating  a  powerful,  miniature  explosion.  Depending  on  details  of  the 
experimental  setup  this  explosion  can  initiate  one  or  more  shocks  in  the 
ambient  gas  as  well  as  initiate  the  turbulent  plume  through  which  the  shock 
propagates  and  interacts  at  a  later  stage. 

Figures  4  and  5  illustrate  the  experimental  pre-shocked  and  post  shocked 
turbulent  power  spectral  density  distributions  and  the  current  LES  simulation 
results.  Amplification  predictions,  particularly  in  the  dominant  low  wave 
number  range,  are  seen  to  be  even  more  satisfactory  than  in  the  previous  low 
Mach  number  results.  The  experimental  decay  is  steeper  than  the  simulated 
decay  in  the  inertial  range  reflecting  an  approach  to  three  dimensionality  (k~3) 
in  spherical  shock  propagation  experiments  in  contrast  to  the  cylindrical 


>* 


Figure  4.  Comparison  of  pre-shocked  normalized  power  spectral 
density  distributions  from  current  LES  computations  with  laser  target 
interaction  hypervelocity  experiments  [1]  at  Mach  numbers  of  the 
order  of  100. 


36 


> 


c  Turbulent  spectra  behind  shock 


Figure  5.  Comparison  of  post-shocked  normalized  power  spectral 
density  distributions  from  current  LES  computations  with  laser  target 
interaction  hypervelocity  experiments  [1]  at  Mach  numbers  of  the 
order  of  100. 


geometry  of  the  simulations  which  are  seen  to  decay  approximately  at  the 
classical  (kt5/3)  inertial  distribution. 

The  shock  amplification  of  the  turbulence  energy  and  the  relaxation 
behind  the  shock  as  well  as  the  distribution  of  the  average  integral  scales  of 
motion  are  displayed  in  Fig.  6,  as  a  function  of  distance  measured  as  fractions 
of  shockwave  radius  of  curvature.  These  distributions  illustrate  perhaps  the 
most  significant  promise  in  the  present  stage  of  LES  development  for 
interaction  analysis.  The  close  parallel  with  the  experimental  results  indicates 
that  crucial  information  can  be  obtained  from  numerical  simulations  on  the 
explicit  grid  scale  range  of  motions,  provided  attention  is  given  to  properly 
modeling  the  influence  of  the  non-resolved  subgrid  scale  motions  on  the 
resolved  motions.  Variation  in  decay  at  substantial  distance  behind  the  shock 
is  the  result  of  physical  differences  between  the  LES  conditions  and  the 
experimental  situation.  The  experimental  post  shock  decay  is  associated  with 
a  rapid  depletion  of  sensible  vapor  matter  near  the  origin.  The  rarefaction  is 
not  simulated  in  the  LES  results  which  were  generated  with  a  constant 
momentum  reservoir  during  the  test  phase. 
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Figure  6.  Comparison  of  turbulent  kinetic  energy  amplification  and 
correlation  length  scale  contraction  through  and  beyond  shock 
interaction  comparing  current  LES  computations  with  hypervelocity 
experiments  [1]. 


4.  SUMMARY 

These  results,  even  though  obtained  at  an  early  stage  in  the  present  trials 
and  comparisons,  seem  to  supply  confirmation  of  the  utility  of  a  formal 
procedure  for  modeling  the  influence  of  the  non-resolvable  subgrid  scales  on 
the  explicitly  computed  scales  of  motion.  At  the  same  time,  however,  recent 
applications  of  adaptive  mesh  refinement  with  the  Godunov  shock  resolution 
schemes  have  produced  much  more  highly  resolved  results  including  a  well 
established  inertial  range  and  appropriate  spectral  decay,  even  in  the  absence 
of  a  sub-grid  scale  model  formal  procedure. 

The  mechanisms  governing  the  influence  of  shock  waves  on  turbulence 
are  not  well  characterized  and,  it  almost  certainly  follows,  not  well 
understood.  Numerical  simulations  and  appropriately  designed  experiments 
may  help  to  reduce  the  uncertainty. 

The  Reynolds  numbers  for  the  particular  flow  circumstances  and 
configurations  of  interest  here  discourage  the  use  of  direct  numerical 
simulations.  Hence,  LES  and  selected  supplemental  computational  techniques 
such  as  the  direct  Monte-Carlo  shock  structure  simulation  procedure  applied 
at  an  earlier  stage  of  this  work  are  the  tools  of  choice  for  augmenting  and 
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analyzing  the  experimental  data.  The  present  effort  focuses  on  current 
progress  in  developing  an  all-Reynolds  number,  compressible  LES  procedure 
that” will  provide  a  valuable  tool  in  investigation  of  the  underlying  physical 
mechanisms  that  govern  the  interaction  of  shock  waves  and  turbulence  in  the 
neighborhood  of  realistically  complicated  geometries. 
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nmnoucTioi 

Tha  interaction  of  two  lnterstreaalng  plaaaaa  la  of  interest  froa 
several  points  of  view.  As  a  basic  pi asaa  physics  topic,  there  are  a 
ayrlad  of  possible  collisionless  and  hydrodynaale  instabilities  that  may 
occur  whan  two  beans  interpenetrate.  Some  aspects  of  these  instabilities 
are  well  established  and  have  been  observed,  but  others  are  poorly 
understood.  Interpenetrating  plasnas  related  to  our  choice  of 
experlaental  paraawters  appear  In  several  natural  and  nan-nade 
situations.  The  sun,  for  instance,  anlta  a  solar  wind  plaaaa  which  foras 
a  flowing  interplanetary  drifting  plaaaa  with  drift  velocities  in  the 
hundreds  of  klloneter  per  second  range;  interestingly,  laser-produced 
plasaaa  have  oonparable  ablation  velocities;  laser-produced-plaaaas  also 
streaa  froa  coronas  having  taaperatirea  vary  close  to  tha  solar  ooronas 
(l-2keV).  The  solar  wind  interactions  with  the  earth's  Magnetosphere  and 
coast ary  pi asass  have  jbeen  observed  to  be  very  coaplex.  An  even  aore 
draaatlo  natural  exaaple  of  interacting  energetic  plaaaaa  la  a  supernova 
explosion.  Han-aade  counter  stressing  plaaaa  situations  also  abound, 
e.g.,  in  soae  proposed  inertial  cooflnaaent  fusion  reactors,  theta  pinch 
devices.  Ion  beaa  I CP  sohsaes,  and  in  soae  ionospheric  Modifications. 

He  describe  a  lasar-produced-plaaaa  experiaent  which  explores 
lnterstreaaing  plaaaa  afreets  in  both  oollislonlesa  and  eollislonal 
reglaes.  One  of  the  plaaaa  ooaponenta  consists  of  the  ablation  plaaaa 
froa  the  Id-laser  solid-target  interaction;  the  other  plaaaa  is  furaed  by 
the  photolonixation  of  tha  aabient  gas  surrounding  the  target.  An 
externally  applied  Magnetic  field  can  be  applied  across  the  interaction 
region. 

He  concentrate  on  the  regiae  where  the  relative  velocity  of  the  two 
plaaaa  coaponenta  have  Mach  niabers  (relative  to  either  sound  or  Alfven 
waves)  auoh  greater  than  one,  and  where  one  coaponent  is  effectively 
unaagnetlzed  while  the  other  coaponent  aay  be  Magnetic  field  doalnated. 
Instability  boundaries,  blast-wave  behavior,  and  Raylelgh-Tsylor  effects 
are  prlae  objectives  of  this  study. 

Hhen  two  plaaaa  coaponenta  interpenetrate  they  streaa  freely  through 
each  other  unless  the  eollislonal  Bean  free  path  is  short  coa pared  to  the 
systea  size.  Collisions! 'coupling  can  occur  through  atoalc.  Molecular,  or 
nuclear  collisions  as  well  as  between  the  various  electron  and  ion 
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.  Collective  plosaj  1  nsiabl  1 1  lies  con  c.tuso  an  "effective* 
llision  frequency  even  when  ordinary  collisions  are  not  dominant.  Beaa- 
i-jij  t nstatt  1  i t ies  can  be  of  the  bump-on-tall  variety,  where  a  small 
'-'am  *  component  sits  on  the  wing  of  the  dominant  component  near  the 
■  ise  velocity  of  plasma  wave;  ours  is  not  this  type.  Our  configuration 
.  of  the  beam-beam  type.  Illustrated  in  figure  I,  where  the  two  counter- 
i -earning  plasma  components  have  a  large  drift  velocity  compared  to  their 
■nTO.il  speeds. 

Strong  interaction  between  the  two  counterstreaming  components  Is 
^served  when  the  colllslonallly  Is  high.  In  this  case,  strong  energy- 
.d-momen turn- conserving  shocks  (also  known  as  ‘blast-waves*  or  Sedov, 

•ylor.  Von  Neumann  shocks)  are  formed.  We  verify  that  the  system  has 
-  .ny  characteristics  predicted  by  a  simple  blast-wave  model,  yet  It  also 
■viates  from  the  simple  model  In  many  respects. 

As  the  ambient  plasma  density  is  reduced,  the  system  becomes 
oillsionless  and  the  coupling  between  the  two  plasma  eounterstreamlng 
components  tend  to  diminish.  However,  signatures  of  collective  beam- 
pi  used  instability,  probably  the  magnetized  ion- ion  (Nil)  Instability, 
are  seen  In  our  experiments.  Although  the  instability  Is  weak  under  our 
conditions.  It  has  the  potential  of  being  an  effective  moment ta  transfer 
mechanism  between  the  energetic  drifting  plasma  and  the  ambient  plasma  la 
i be  colllslonless  regime. 

We  describe  the  behavior  of  the  interaction  found  experimentally  and 
rake  comparisons  where  possible  with  theory  and  hydrodynamic  code 
alculatlons. 

nXPERIHENT 

The  experiment  consists  of  focusing  the  beam  from  the  HRL-Pharos  II 
Kd-laser  (1.05  pm  wavelength)  onto  small  foil  targets  (£  1  am  dla,  few- 
ml crons  thick  A1  or  CH)  In  the  center  of  the  target  chamber.  Typically, 
the  laser  pulse  Is  a  gausslan  *-na  FHfll  100  J  pulse  focusad  In  the  range 
of  1012  to  10'*  W/cm2.  Target  material  is  ablated  by  the  laser 
irradiation  and  streams  radially  outward  at  hlgh-vmloclty  (20*- 100* 
km/s).  We  call  this  plasma  component  the  target-plasma.  A  low  density 
background  gas  maintained  the  chamber  is  promptly  ionized  In  the  vicinity 
of  the  target  by  radiation  from  the  laser-target  interaction 
(photolonixatlon),  and  additionally,  at  later  tines  by  the  expanding 
target-plasma  (UV  photolonixatlon  from  the  plasma  shell  emission  or  by 
oartlcle  Impact).  The  background  gaa  is  usually  nitrogen  although  gases 
from  the  hydrogen  to  Xm  have  also  been  used.  This  creates  a  stationary 
ambient  plasma  through  which  the  high- velocity  target-plasma  streamsi  an 
external  magnetic  field  Is  sometimes  applied  over  the  whole  interaction 
-eg Ion.  using  permanent  magnets  or  helaholts  configuration  oolls.  The 
spatial  scale  of  the  Interaction  region  Is  centimeters,  the  time-scales 
involved  are  tens-of -nanoseconds,  the  magnetic  fields,  when  applied,  are 
in  the  kllogauss  range,  and  the  ambient  gas  pressures  used  are  ten  Terr  (1 
Torr  -  1  an  of  Hg  -  3.2*  a  10* 6  aolec/cc)  and  below.  Table  I  summarizes 
the  experimental  parameters  and  Figure  2  shows  the  experimental 
arrangement  schematically. 

DIAWOSTICS 

Many  Instruments  are  used  to  measure  the  evolution  of  the  beam-plasma 
Interaction.  Diagnostics  include:  Incident  and  reflected  laser-beam 
energy  calorimeters  and  Use- re  solving  light-diodes  to  measure  the  time- 
history  or  the  Incident  laser  pulse.  The  angular  distributions  of  the 
resulting  target-plasma  energy,  velocity,  mass,  and  momentum  are  obtained 
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Figure  1.  Schematic  of  the  target-planaa  shall  Boring  outward  froa  th« 
laaar  focal  region  though  a  magnetized  photo Ionized  aablent 
gas  (left).  The  aablent  (stationary)  plasaa  and  drifting 
target-plaaaa  (ablation)  distributions  fora  s  classic 
lnterslreaalng  instability  configuration  (right). 
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Parameters  used  in  laser -plasma  counters! reaming  experiment 


Laser : 

Energy 

Pulse  Duration 
I rr ad lance 

Ambient  Gas: 

Pressure 

Species 

Target  Plasma; 

Initial  Velocity 
Hass  (At) 

Magnetic  Field: 


l-aoo  J 


«-ns  FWHM 
10'3- 


Torr 

H,  He.  N2,  At.  He.  Air,  Xe.  Kr... 


100-1000  Ka/ sec 

0.1-1. Out* 

0-1000  gauss 


with  arrays  of  alnl-calor testers,  tlae-of-fllght  Ion  detectors,  and  light 
sensitive  diodes.2  Magnetic  loop  Induction  probes  aeasure  Magnetic-field 
distortions  in  the  target-plasaa/aablent-plasaa  interaction  region  and  any 
self-generated  magnetic  field  components  as  a  functions  of  position  and 
ilee.’  Optical  laser-probe  diagnostics,  such  as  dark-field,  two-tlae 
shadowgraphy  and  lnterferoaetry,  yield  quantitative  pictures  of  the  plasaa 
density  structure  at  aoderate-to-hlgh  densities  (>  10'7  electrons  per 
cn3).  Laser  scattering  provides  density  and  teaperature  Information  froa 
Thoason  scattering,  as  well  as  plasaa  fluctuation  spectral  Information  In 
tne  collective  reglrs.  Optical  laaglng  with  still  photography,  time- 
resolved  framing  photography,  and  streak  photography  also  give  qualitative 
and  quantitative  Information  on  plasaa  properties  and  structure. 
Spectroscopy,  froa  the  Infrared  to  XUV,  is  used  to  give  spatially  and 
temporally  resolved  Information  about  the  velocity,  temperature, 
ionization  state,  density,  emlsslvlty,  opacity,  etc.  of  the  target-plasma 
material,  the  photolonlsed  background  plasma,  and  the  Interaction  between 
them.’  X-ray  and  XUV  diagnostics  Alow  diagnosis  of  plasma  properties 
above  20  eKj  X-ray  pinhole  photography  gives  two-dimensional  apatlally- 
resolved  Images  while  electronic  s-ray  and  UV  detectors  give  quantitative 
information  about  the  Initial  plasma  temperature,  the  radiant  emission 
from  the  expanding  tar  get- plasma/ ambient- plasma  rront.  Tracer-dot 
techniques.  In  which  small  spots  of  a  hlgh-Z  aaterlal  are  Implanted  In  the 
target  surface,  have  proved  valuable  In  providing  target -plasma  flow 
visualizations  and  spatially  resolved  quantitative  spectroscopy.6 

PLASMA  PROPERTIES 

The  laser-produced  target-plasma  conditions  are  characteristic  of  the 
steady-state  laser-plasma  ablation  process. 2,7  This  plasma  Is  well  suited 
for  this  experiment  since  It  has  a  reproducible  single-peaked  velocity 
distribution.  The  mean  velocities  of  the  Ion  distributions  Vd,  are 
tunable6  simply  by  varying  the  laser-irradlance  Xi  Vd  varies  with  the  0.2 
sower  of  I.  A  good  example  of  an  Ion  tlae-of-fllght  trace,  from  which  the 
velocity  distribution  of  the  target  plasaa  can  be  determined,  is  shown  In 
Figure  3>  In  addition,  the  velocity  spread  is  relatively  narrow, 
typically  AV/V  Is  about  0.25;  AV«  >  is  also  tisiable  from  0.1  to  1  by 
increasing  the  focal  spot  disaster  to  pulse  duration  ratio.9  The  mass 
ablation  rate  of  the  target  plasaa  Is  also  well-known,  and  goes  as  the  0.6 
power  of  laser  Irradlance. 
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FI pjr*  3.  Laser  target-plasma  Ion  tl*e-of-fllght  detector  signal  (left) 
and  velocity,  moment;*  and  energy  distributions  (right).  Note 
the  high-velocity,  narrow-velocity  spread  target-plasma. 

These  ions  contain  over  60S  of  the  absorbed  laser  energy. 


The  photolqnlxed  aablent  plasma  properties  have  been  primarily 
determined  by  using  plasma  spectroscopic  methods.  These  techniques  have 
been  described  elsewhere.’1 1 0  Typical  ambient  plasma  temperatures  are  in 
the  1.5  to  2.5  eV  range,  depending  upon  the  gas  density.  The  degree  of 
ionisation  is  high  near  the  laser-target  focal  region  and  decreases  to  a 
few  tenths  of  a  percent  about  1-em  away.  Of  course,  when  the  target 
plasma  eventually  streams  through,  the  ambient  properties  change 
dramatically;  these  changes  will  be  described  in  a  later  section. 

COLLISIONS!.  INTERACTIONS 

A  strong  shock  rarms  when  an  energetic  plasma  burst  expands 
supersonically  into  another  plasma  when  the  collision  me  an- free- path  is 
small.  The  shock  propagates  Into  tha  ambient  plasms,  sweeping  it  up  into 
a  thin  coupling  shell,  which  consequently  slows  down  due  to  the  mass 
accretion.  If  the"  Initial  energy  Is  released  quickly  compared  to  the  'lime 
scales  of  Interest  and  both  particle  energy  and  momentum  are  conserved, 
the  resulting  shook  front  Is  termed  a  Taylor-von-Neumann-Sedov  shock  or  a 
•blast-wave.*1 1 

In  this  section  we  review  some  features  or  the  blast-wave  model  and 
use  than  to  Interpret  the  properties  of  coupling  rronts  observed  In  the 
experiment.  He  find  good  agreement  between  experimental  results  and  most 
blast-wave  theory.12  However,  in  contrast  to  an  Ideal  blast-wave,  which 
Is  hydrodynamlcally  stable,  the  shocks  In  the  laser-experiment  develop 
striking  spatial  structure,  resembling  arterial  aneurisms,  under  certain 
circumstances.  Causes  of  these  nonunirormltles  are  not  yet  Isolated; 
none-the-less,  we  speculate  on  some  possible  responsible  mechanisms. 

The  parameters  used  in  the  eolllslonal  regime  experiments  are 
Included  In  Table  I.  Notice  that  the  parameters,  were  varied  over  a  broad 
range  to  adequately  test  the  blast-wave  aodel  scaling.  Also,  In  some 
shots  a  500  C  magnetic  field  was  applied  over  the  interaction  volume 
(transverse  to  the  laser  beam).  However,  no  magnetic  field  dependence  was 
seen  In  the  eolllslonal  regime.  Dark-field  shadowgrams  were  taken  of  the 
shock  structure  at  several  times  after  the  laser  pulse.  Spectroscopic 
observations  were  also  made  to  determine  the  state  (density  and 
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*»r  i  h*»  tmbifinl  .iful  i.’oupled  plasmas. 


.'tamper  oi  jl.1^  has  seen  numerous  examples  or  coupling  shells  (shock 
ml),  such  as  the  one  shown  In  Figure  k,  taken  with  dual-liae  dark-field 
■er-probe  shadowgraphy.  These  photographs  indicate  that  the  shells  have 
.<■  following  general  features: 

*  A  thin  (AR/R  *  0.03),  approximately  spherical  shock  is  observed 
propagating  into  the  aablent  aedia  at  tlaes  long  after  the  laser 
pulse  has  teralnated. 

*  The  shocks  decelerate  as  they  propagate  away  froa  the  focal 
region. 

*  The  velocity  of  the  shock  is  a  function  of  the  deposited  laser 
energy,  aabient  gas  type  and,  of  course,  tlaei  but  the  Motion  or 
the  shell  is  Insensitive  to  the  initial  target-plasaa  velocity. 

However,  the  shocks  also  develop  structure,  such  as  shown  in  Figure 
5.  at  the  higher  aabient  pressures.  We  shall  return  to  this  point 
later.  The  unperturbed  portions  of  these  shells  follow  the  aaaM  blast- 
wave  scaling  as  totally  unperturbed  shells,  but  the  spatial  perturbations 
deviate  froa  the  blast-wave  behavior. 

The  spectroscopic  results  indicate  that  the  aabient  plaaaa  is 
initially  weakly  ionised  (0.2f)  at  1-2  eV,  one- centimeter  froa  the  target 
surface.  But  when  the  blast-front  arrives,  the  plaaaa  becomes  1001 
ionized  with  a  temperature  of  about  10  eV;  the  mass  density  Jiasps  above 
me  initial  aabient  level  by  a  factor  of  7-10  in  the  shock. 


Figure  A.  Dual-tlae  dark-field  shadowgraa  of  shook  waves  at  52  and  96 
nsec  In  5  Torr  ¥»  gss.  The  incident  laser  energy  was  A.1  J 
and  the  initial  debris  valoclty  was  2  s  10T  ca/see;  B-0. 
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figure  5.  Dual-time  dark-field  shadowgrams  of  shockwaves  with 

aneurisms.  (tert)  Shadowgram  of  a  shock  wav*  at  52  and  96 
naac  In  a  5  Terr  anbiant  (90S  M2  *  10%  H-)  gas.  The  laser 
energy  was  38  J  and  the  initial  target-plasma  speed  was  5  x 
10'  cm/seci  8*0.  Note  the  growing  "aneurism*  at  the  9:00  pa 
poelllon.  The  object  on  the  right  Is  a  nagnetlc  probe  (out  of 
focus).  (Sight)  Shadow gran  of  shock  fronts  in  a  1.5  Torr 
(Mj  ♦  Ho)  gas.  The  observation  tlaes  were  at  52  and  16*  nsec, 
the  Incident  laser  energy  was  20  J,  the  initial  debris  speed 
was  approxlaately  3  x  107  ca/sec,  and  a  600  gauss  nagnetlc 
field  was  present  Into  the  plane  of  the  paper.  The  gaps  In 
the  target  holders  are  about  5  mm. 


Ve  shall  coapare  these  experimental  observations  with  a  blast-wave 
■odel.  After  the  Initial  energy  release  (laser-pulse),  the  plasma  rapidly 
expands,  picking  up  ambient  eedla  along  the  way.  After  the  shell  has* 
accreted  an  —blent  mass  several  tines  the  Initial  target-plaaaa  mass,  the 
shell  decelerates  with  the  familiar  self -stellar  blast-wave  dependence 
*«(E/p),'5tz  5.  Eventually,  when  the  shell  velocity  approaches  ths 
acoustic  speed  In  thf  ambient  media,  the  disturbance  is  no  longer  shock- 


treatments  of  blast-waves  since  Taylor,  Von  Neumann,  and  Sedov.  Some  of 
these  works  extend  the  theory  into  the  initial  phase,  where  the  target- 
plasma  mass  Is  Important,  while  others  are  hydrodynamic  ealculatlisa. 

Here,  we  follow  the  method  or  Chernyl  as  outlined  In  Zaldovlch  and 
Nalzer.  This  blast-wave  approximation  has  been  shown  to  yield  results 
within  a  few  percent  of  exact  treatments.  The  following  assumptions  are 
made i  1.  The  energy  release  Is  considered  an  instantaneous  point 
explosion.  2.  Spherical  symmetry  Is  assumed  for  simplicity.  3.  The 
laser-target  velocity  and  the  resulting  shock  speeds  are  much  larger  than 
the  undisturbed  ambient  sound  speed.  *.  The  expansion  conserves  particle 
energy  and  momentum.  5.  The  ambient  gas/plasma  Is  swept  up  by  the  debris 
front  Into  a  thin  cold  shell  having  a  mass  large  compared  to  that  of  the 
Initial  laser-target  plasaa.  6.  The  media  is  characterized  by  a  constant 
effective  ratio  or  specific  heats  V.  7.  Finally,  counter-pressure  due  to 
the  ambient  plasma  Is  neglected. 


The  shell  front  Is  a  strong  shock  wave  and  the  Hugonlot  jimp 
relations  apply  between  the  ambient  media  (6)  and  shell  (s).  The  density 
Jimp  Is  therefore  given  by 
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T*ie  flow  velocity  behind  the  shock,  u9.  Is  related  to  the  shock  speed  V# 

'•y 


riving  the  pressure  within  the  shock. 


%  '  <£>  P0 Vf- 


(2) 


(3) 


Now.  combining  these  results  of  the  strong  shock  Jump  relations  with 
conservation  of  sass,  energy  and  momentim  we  obtain  aany  blast-wave 
woperties.  Conservation  of  aass  is  expressed  by. 


k*R24Rp. 


U»RJ 


p.  ♦  («,)  -  H. 


{*> 


The  quantity  on  the  left  side  of  Eqn.  (A)  Is  the  total  shell  lusa  as  a 
function  of  shell  radius  R  and  shell  thickness  ARs  on  the  right  side  Is 
the  aass  of  ambient  gas  within  the  bubble  voluae  ( assumed  to  be 
completely  swept  up)  plus  the  initial  debris  aass  a^  (neglected  here). 
The  relative  thickness  of  the  shell  Is  found  by  combining  Eqn.  (A)  with 
Eqn.  (1),  i.e.. 


AR  1  T — 1 
R  '  3  Y«1 


(5) 


Proceeding  further,  conservation  of  eoaentua  is  expressed  by, 

&<"“.*>  ■  (6) 

?b  is  the  pressure  within  the  bubble  voluae  which  pushes  outward  on  the 
shell.  The  shell  Is  assiawd  to  have  aost  or  the  system  aass,  but  scan 
mall  amount  of  aass  must  remain  Inside  the  shell  boundary  (bubble). 
~inally,  conservation  of  energy  sets  the  energy  In  the  explosion  E  equal 
to  the  sum  of  the  shell  kinetic  energy  plus  the  thermal  energy  Invested  In 
the  system;  this  Is  expressed  asi 

E  ‘  \  Muf  *  t=7  ^  pb  *  (t=7  *»»2«V*  ») 


The  first  term  on  the  right-hand-side  of  Eqn.  (7)  Is.  the  shell  kinetic 
energy,  the  second  and  third  teras  are  the  theraal  energies  within  the 
bubble  and  shell  respectl vely.  The  last  term  (shell  theraal  energy)  Is 
usually  neglected  relative  to  the  second  tera  (bubble  theraal  energy) 
since  the  ratio  is  or  order  10  .  -  Me  aake  the  same  assuaption  here, 
however,  note  that  these  two  contributions  to  the  theraal  energy  beccae 
sore  comparable  as  AR  Increases,  as  Pb  decreases,  or  In  the  event  that 
the  v  of  the  plasma  in  the  bubble  Is  higher  than  that  of  the  shell  (which 
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could  be  true  since  the  bubble  has  a  such  hotter,  lower-density  plasma 
than  the  shell).  This  assumption  about  the  apportionment  of  thermal 
energy  does  not  change  the  blast-wave  radlis  versus  time  scaling  -  only 
the  constant  of  proportionality.  Ue  shall  return  to  this  point  again. 
Demanding  that  the  energy  E  be  Independent  of  radius  throughout  the 
expansion  and  assuming  that  P„  «  P#  gives  a  bubble  pressure  about  half 
that  of  the  shell  pressure  P.,  i.e. , 


P  *  —  P 

rb  2  V 


(8) 


[This  is  to  be  compared  to  P#  •  O.li  P,  for  T  •  1.2  in  the  exact  case.] 

How,  from  the  above  relations,  the  expression  for  the  blast-wave 
radius  with  time  similarity  solution  is 


*(E.p0.t) 


«o  <EV,/5t2/5, 


(9a) 


or,  in  "practical*  units. 


R(cm)«0.092  coCE{J)/(P(Torr)/(«f/MW|(  )) J1/5t(nsec)2/5.  (9b> 


where  )  is  the  ratio  of  the  ambient  gas  molecular  weight  relative 

to  a  nitrogen  ^molecule,  and  (0  is  a  weak  function  of  Y  with  a  value  of 

order  unity.  Hlthln  our  met  of  assiaiptlona,  c  is  given  by  the  relation, 

o 


75  (T-1)(Y«1)  .1/5 


co  ■  (i6,  (iv-n — * 


(10) 


For  completeness,  we  extend  the  treatment  In  Ref  (11)  to  Include  the 
shell  thermal  energy  In  the  energy  balance  [third  term  of  Eqn.  (7)).  2  Me 
also  allow  far  the  Y  of  the  plasma  within  the  bubble  to  differ  from  the 
shell/amblent  plasma  Y  by  designating  the  bubble  y  by  Y.  and  that  of  the  > 
remaining  plasma  by  Y  i  then  Eqn.  (10)  becomes. 


(Y„-1)(Y"1)2 
o  _  ■% 


1/5 

m 


(10‘) 


The  ratio  or  inferred  explosive  energy  release  under  the  two  sets  or 
assumptions  (g  /r  1 y  can  differ  by  about  a  factor  of  two  although  the 
maximum  error  ?n  8(t)  Is  only  125.  It  is  clear  that  detailed  hydrodynamic 
calculations  which  keep  track  or  the  local  values  of  Y  are  necessary  to 
get  a  precise  description  of  the  expansion.  He  use  Eqn.  (10)  in  the 
remainder  of  this  paper. 

The  ratio  or  thermal  energy  to  kinetic  energy  in  the  blast-wave 
system  Is  surprisingly  high;  this  ratio,  obtained  by  taking  the  ratio  or 
the  second-term  to  first-term  In  Eqn.  (7),  is  given  approximately  by 


«EE 


.  1  (El). 

2  Vr-r* 


(ii) 
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[  Tne  right  hand  side  of  Eqn.  (11)  becomes  ( i /?)[  (v*1  J/df^-l )  1 
■inder  1  he  same  sel  of  assumptions  as  Eqn .  (10').]  Other  relevant  blast- 
wave  parameters  are  the  plasma  effective  v  and  the  temperature  In  the 
snell  and  in  the  bubble  volume.  The  temperature  In  the  shell  can  be 
estimated  by  using  an  approximation  to  c,  the  internal  energy  of  air. 


t  -  8.3  Ta(eV)'*5  (pA/ps)0*12  eV/molec, 


(12) 


which  is  valid  for  temperature  T#  between  1  and  25  eV,  and 
density  p  between  <0p.(p.  •  atmospheric  density)  and  l0“3/p.st  ranges 
from  1.1  fo  1.3  for  air  in  this  regime  with  Y  *  1.2b  a  good  •effective'* 
value.  The  internal  energy  is  also  given  by 


t 


1  P 
T-1  »' 


(13a) 


where  P  and  p  can  be  determined  through  Eqns.  (1),  (3).  (6)  or  by  direct 
measurement.  In  the  shock  front  Eqn.  (13a)  becomes 


c  - — - — y  Joules  per  kg,  (13b) 

(T*1)Z 


or,  to  obtain  Eqn.  (13b)  In  the  sane  units  as  Eqn.  112),  multiply  by  0.33b 
x  HU  and  express  the  shock  speed  Ya  in  units  of  (107  cm/ sec).  Equating 
Eqn.  (13a)  to  (12),  with  appropriate  units,  gives  an  estimate  for  T_.  The 
resulting  expression  for  temperature  in  the  shell  is  thereby  found  to 
be. 


T#(eV) 


v^(.i07  cm/sec) *hu1  2/3 

.  <T«1)2(pA/ps)0,12J 


(lb) 


A  tabulation  of  these  blast-wave  parameters  is  given  in  Table  II 
for  Y  •  1.2,  l.b,  and  5/3.  Dote  that,  as  assuaed,  most  or  the  mass  is  in 
a  very  thin  shell.  Also,  the  high  temperature  within  the  bubble  is  a 
consequence  of  the  approximate  pressure  balance  with  the  shell  (but  with  a 
much  lower  density).  As  we  go  towards  the  center  or  the^bubble. .the  model 
predicts  that  the  plasma  density  goes  to  sero  ast  p  -  R3  t  , 

a  *  f»  V  • 


Blast-wave  Experiments 

We  now  compare  the  experimental  findings  with  the  blast-wave  model. 
The  main  observables  in  this  experimental  series,  that  we  will  relate  to 
blast-wave  theory  are  the  shell  position  R,  the  thickness  of  the  shell  AR, 
and  density  p  and  temperature  T#  of  the  shock  wave.  Experimental 
variables  included:  the  laser  energy,  the  laser  focal  spot  sice  (and 
thereby  the  initial  target-plasma  velocity),  the  ambient  gas  type  and 
pressure,  the  presence  or  absence  of  an  external  magnetic  field, 
occasional  variations  in  the  target  angle  or  structure,  and  the 
observation  times. 

Shell  position  and  blast-wave  scaling.  A  plot  of  the  distance  of  the 
shock  fronts  rrom  the  target  surrace,  R,  For  experimental  shots  which  span 
the  range  of  parameters  tabulated  in  Table  I,  is  shown  in  Figure  61  the 
variables  along  the  abscissa  of  Figure  6  are  scaled  according  to  Eqn. 
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Figure  6.  Plot  of  shock  from,  positions  R  as  s  fucntlon  of  the 

normal 1 ted  blast-wave  scaling  parameter  for  the  data  set  in 
!•  note  the  excellent  consistency  with  blast-wave 
scaling  with  0.092  «0  -  0.123. 


Note  ihe  good  agreement, of  the  enlire  data  sel  with  the  blast -wave 
scaling  parameter  [(E/p  H‘j  .  with  a  single  universal  constant  of 
proportionality,  t  -  8.123/0.09?  •  1.3*  (from  Eqn.  9b).  2  The  scaling 
is  insensitive  to  Pfie  Initial  target-plasaa  velocity  for  constant  incident 
laser  energy. 


It  is  possible  to  relate  the  experimentally 
observed  (  to  theory  (  by  taking  into  account  the  angular  distributions 
of  the  plaSSa  expansion.0  If  we  use  the  fact  that  half  the  laser-plasaa 
energy  is  contained  within  a  half-cone  angle  or  *0*  from  the  normal  of  the 
target  in  vacuum  and  assume  that  this  angular  distribution  still  holds 
true  thoughout  the  expansion  (thia  may  not  be  too  bad  an  assumption  since 
the  flow  is  very  supersonic),  then  we  obtain  an  equivalent  spherical 
experimental  value  for  of  1.0  *  0.1. 

Coupling  Efficiency.  Ho  distinct  laser-target-plaaaa  ion  peak 
reaches  our  tlme-of-fllght  detectors  at  2  Torr  fill  pressure,  and  most  of 
the  target-plasaa  ion  peak  is  lost  at  200  aTorr.  He  conclude,  t  her  fore, 
that  the  coupling  between  target-plasma  and  ambient  plasma  is  high  in  this 
pressure  regime,  and  nearly  complete  above  1  Torr. 


The  coupling  efficiency, 
dependence  can  be  inferred  by 


^lajt^w?ye/Elaser  • 


ave- -laser*  »“61*nt  PrMaur« 
the  radlus-tlae  observations  to  the 


blast-wave  model  and  solving  Tor  E-E^  in  Eqn.  (9).  Figure  7  shows  the 
results  for  an  initial  target-plasaa  velocity  of  A. 5  x  10'cm/s  streaming 
into  nitrogen  ambient  plasma.  The  coupling  efficiency  drop  below  1  Torr 
is  roughly  consistent  with  the  loss  of  colllstonallty. 


Shock  Thickness.  The  shell  thlekness-to-radlus  ratio  AR/R  is 
observed  to  be  about  0.03  l  0.01.  In  fact,  the  brlght-dark-brlght 
structure  seen  in  the  shock  front  shadowgrams  are  indicative  of  a  steep 
gradient  on  both  the  front  and  back  surfaces  of  the  shell.  This 
implies  T  ■  1.20  t  0.07  from  Eqn.  (S).  a  value  consistent  with  both  the 
determination  from  R(t)  (above)  and  t  T>  tion-of-state  or  air, 

Eqn. (12).  Actually,  the  shell  thickn  a  relatively  sensitive 

independent  indicator  of  the  effective  Y;  for  convenience,  wo  invert  Eqn. 
(5)  and  solve  for  Y,  l.O. , 


Y 


1*3(AR/R) 

TTramn- 


(15) 


McLean  et  ml.5  uses  results  from  spectroscopic  continuum  measurements 
to  infer  the  density  of  the  plasma  within  the  shock  front.  Typical  shell 
densities  are  found  to  be  about  7  to  10  times  the  ambient  N2  density  above 
1  Torr  fill  pressure.  The  Inferred  Y  from  Eqn.  (1),  expressed  by 


a/>o) 

■(V^r=r* 


yields  Y  •  1.18  t  0.04  using  the  observed  density  Jumps.  Thus,  the 
density  Jump  at  the  shell  la  also  consistent  with  a  blast-wave 
with  Y  -  1.2  and  the  other  experimental  results. 

Interferograas  of  the  shock  confirm  the  shock  density  Jump  obtained 
spectroscopically.  Figure  8  shows  such  a  shock  profile  obtained  from 
interferometry.  The  interferograas  also  show  behavior  not  predicted  by 
the  simple  blast-wave  model.  Hon-blast-vave  features  observed  include;  a 
density  reap  and  step  in  front  of  the  steep  shock  bubble  and  higher  than 
expected  density  Inside  the  shock  bubble.’3  Moreover,  aneurisms  are 
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Figure  7.  Coupling  efficiency  of  the  interaction  between  a  t.5  x  10‘ 
ca/s  laser-target  plasaa  and  aablent  nitrogen  plaaaa. 


Figure  8.  Electron-density  shock  profile  obtained  froe  an 

Interferograe.  Mon-blast-wave  like  features  are  noted. 
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-l»F  mitely  not  m  tnc  blast-wave  loouel .  Figure  9  shows  another  dual-tine 
■.n.ioowgram  showing  Further  non- blest-wave  Features.  In  addition  to  the 
.seueisffl,  considerable  plasaa  turbulence  or  FluFF.  Is  seen  to  the  rear  oF 
me  target  (away  froa  the  laser  beaa).  This  structure  is  probably  caused 
:>y  the  disruption  or  the  target  aaterlal  accelerated  by  the  ablatlon- 
;>lanoa.  Note  also  in  Fig.  9  that  there  is  scattered  laser-probe  light  In 
me  bubble,  on  the  inside  oF  the  blast-Fronls;  this  is  likely  to  be  caused 
■•y  short-wavelength  plasaa  turbulence  in  the  bubble. 

Shock  Teaoerqtures.  Shock  front  teaperatures  of  10-15  eV  were 
estlaated  by  McLean  et  al.“  froa  the  highest  ionisation  state  or  nitrogen 
observed.  The  teaperature  predicted  within  the  shock  front 
[using  y  •  1.2,  p./p  *  100,  MW  •  26,  and  typical  shock  speed  at  6-1  ca  of 
v3  •  1  (x  107  ca/8es7  in  Eqn.  (16)]  is  Ta  •  9  eV.  This  is  also  in 
reaarkable  agreeaent  with  experiaent. 

The  very  low  density  plasaa  within  the  shell  cavity,  or  bubble, 
should  be  at  a  auch  higher  teaperature  than  T#.  Ho  aeasureaent  of  Tb  has 
yet  been  aade.  But,  to  estlaate  what  to  expect,  we  asauaa  that  the 
equatlon-of-state  of  this  plasaa  continues  to  follow  Eqn.  (12)  [not  too 
likely  since  Eqn.  (12)  is  baaed  upon  Saha  equilibria  and  the  bubble  should 
be  closer  to  coronal  equllibriua],  then  the  bubble  teaperature.wlll  be 
higher  than  the  shell  teaperature  by  a  factor  of  order  (p  /pb)u"  * 
scaling  was  obtained  by  ass  using  pressure  balance  throughout°the  blast- 
wave  systea,  which  sets  cp  -  constant.  A  acre  accurate  air /tar get -plasaa 
equatlon-of-state  for  the  bubble  plasaa  is  needed  to  Bake  a  better 
estlaate.  Measureaent  of  Tb  is  an  axperlaental  challenge  due  to  the  low- 
density  of  the  bubble  plasaa  within  the  hlgh-dcnslty  shell. 


Figure  9.  Dual-tine  shock-front  shadowgraas  (55  and  160  ns)  of  a  5T  H, 
sab lent  shot  with  36  J  on  a  stalk  aounted  A1  target.  Note  tl 
blast-wave  fronts  and  the  non-blast-wave  like  turbulence  and 
aneurisa. 
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Shock  From  jjomwlfOTillM.  Whal  oajses  ihc  shock  from 
nonun 1 focal ties  that  art  observed  to  develop  in  the  exper  nent?  Why  .v* 
the  "aneurism*  nonuni  formulas,  such  as  seen  in  Figures  S  and  9,  so 
uierd?  To  answer  these  questions  we  need  Inventive  theory  and  experiments 
to  ellainatt  or  conflrs  aechanisas.  Here  are  a  few  speculations. 

Expanding  ideal  shock  waves  tend  to  he  hydrodynaalcally  stable,  yet 
this  statement  has  not,  to  our  knowledge,  been  proven  in  general.  IT  the 
shock  fronts  are  Raylelgh-Taylor  unstable  for  some  reason,  (caused  by.  Tor 
exaaple,  an  adverse  density  gradient  set  up  by  radiation  energy  loss  froa 
within  the  swept  up  debrls/aabient  plasaa)  the  growth-rates  can  be  very 
large.  For  exaaple,  taking  shell  decelerations  typical  of  the  experiaent 
(g  -  5  x  101*  <m/sec?)  and  typical  wavelengths  observed  (1-3  an)  yields 
growth  rates,  T__*  (kg) 1/2  of  order  io8/sec,  sufficient  to  create  large 
non  uniformities  within  typical  expansion  Uses. 

Another  possible  mechanise,  target  Jetting,  could  cause  aneurism-type 
protrusions.  Buspe  in  the  coupling  front  night  occur  due  to  the  impact 
of  slower  target  debris  with  the  decelerated  blast-wave.  But,  aneurisms 
still  showed  up  in  an  experimental  aeries  using  thin  foils  and  llaited 
aass  targets  which  should  have  been  completely  ablated  by  the  laser  pulse, 
thereby  eliminating  a  source  of  slower  debris  material.  The  aneurisms 
often  occur  outside  the  incident  laser- be  am  path  so  that  incident  laser 
beam  effects  are  excluded.  Keskinen5  has  proposed  interesting 
asyanatrlxing  aechanisas  caused  by  the  self-generated  magnetic  fields  that 
aay  be  preeent  during  the  initial  expansion;  these  magnetic  fields  modify 
the  flow  patterns  of  the  expanding  debris  plasms.  A  novel  explanation  by 
J.  Giuliani, 6  also  saploys  the  self-generated  Magnetic  fields  trapped 
inside  of  the  shell.  In  this  model  hot  plsaaa  Inside  the  bubble 
prefer entially  sblates  the  back  of  the  central  region  of  the  blast-front 
near  tbs  axis  of  symmetry  due  to  Magnetic  Insulation  off  axis,  thereby 
pushing  out  an  aneuriaa  near  the  target  normal. 

Other  nonuni foralty-induclng  aechanisas  arc,  no  doubt,  possible;  a 
full  understanding  or  this  phenomena  awaits  further  study.  It  Is  noted 
that  some  similar  shock  front  nonuniroraltiss  have  also  been  seen  by 
others.1' 

Ms  have  seen  that  highly-coupled  blast-waves  are  formed  at  pressures 
above  0.$  To rr  In  the  laser- experiment.  These  shocks  are  thin 
(AR/I  9  0.03)  dens#  (p  / p  9  10)  cool  (T  9  10  eV)  and  exhibit  many 
properties  associated  wltn  energy-  and  moment  urn- conserving  blast-wave. 
However,  considerable  non-ideal  non-blast-ware  features  appear  upon  dose 
Inspection. 

C0U4SX0NLESS  'INTERACTIONS 

It  Is  evident,  from  Figure  7,  that  the  Interaction  between  the  fast 
laser-produced  target-plasms  and  the  stationary  ambient  plasma  diminishes 
as  the  ambient  pressure  la  reduced.  This  is  primarily  due  to  classical 
collision  aeaa-fres-pmtha  becoming  large  compared  to  the  Interaction 
region  and  tbs  two  plasmas  cease  Interacting.  In  this  colllaionless 
regime,  however,  plasms  instabilities  can  occur  and  still  causa  momentum 
exchange  between  the  two  plasaa  eoaponenta;  this  gives  an  effectively 
higher,  or  ’anomalous",  collision  frequency. 

In  our  first  experiments  in  the  low-pressure  (<200  nTorr)  regime,  we 
searched  Tor  signs  or  beam-plasma  instabilities.  Five  streaming 
instabilities  are  thought  to  be  potentially  important  In  our  experimental 
regime;18  these  Instabilities  are:  the  magnetised  and  unaagnetlxad  ion- 
ion.  modified  two-stream.  Ion  sound,  and  beam  cyclotron  Instabilities.  Ve 
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!•»■•«»  our  alloni ion  on  the  magnetized  ion-ion  (HID  instability.  The 
"Utfncl i zou  ion-ion  instability  is  an  officiant  momentum  transfer  mechanism 
aitn  jn  effective  collision  frequency  of  about  one-tenth  the  lower  hybrid 
frequency. 


He  are  in  a  regime  in  which  the  target  plasms  Ions  are  magnetized 
-cyclotron  radius  >>  equal  aass  radius),  the  aablent  plasaa  ia  magnetized 
> electron  and  ton  cyclotron  radii  <  the  equal  aass  radius),  and,  aosl 
importantly,  the  target-plasm*  aagnetlc  Mach  nuaber  Is  high  (V.  > 

V»lfven)- 

The  uncoupled  target-plasma  ion  distribution  has  high-velocity  (up  to 
8  x  10'  ca/sec)  and  a  narrow  speed  as  noted  previously.  Thus,  the 
relative  velocity  between  target-plasaa  and  aablent  ion  distributions 
exceeds  the  thermal  speeds  of  the  beams,  satisfying  one  of  the  criteria 
for  streaming  instability. 

Healdual  Collisions.  When  no  magnetic  field  la  applied-,  but  ambient 
gas  is  present  below  80alllltorr,  the  target-plasma  Ion  distributions  are 
not  generally  altered.  For  ambient  pressures  of  80  nlllltorr  and  up,  the 
target-plasaa  Ion  distributions  reaching  the  tlae-of-f light  detector 
(about  22  cm  away  from  the  target)  are  somewhat  attenuated  and 
broadened.  These  pressure  effects  on  the  target-plasma  velocity  Ions 
distribution  are  thought  to  be  due  to  colllslonal  processes,  such  as 
charge  exchange,  since  they  do  not  seem  sensitive  to  the  presence  or 
absence  of  a  aagnetlc  field  (at  least  not  sensitive  to  fields  below  one- 
kllogauss). 

This  hypothesis  was  confirmed  by  using  three  time-of-flight  charge 
collector  detectors  placed  10,  25,  and  55  centimeters  from  the  target,  to 
detect  the  ion  distribution  change*  occurring  between  them. 

One  question  is  whether  we  can  rely  upon  the  preservation  of  the 
shape  of  the. ion  distribution  in  passing  through  the  aablent  media?  If  so 
we  can  detect  the  effects  of  an  interaction  occurring  within  a  few 
centimeters  of  the  target  much  further  away.  The  answer  is  yes,  at  least 
through  masa-pathlengtb  of  5  x  I0,b  molecules/ co  •  os  (l.e.,  150  aforr  », 
with  t  •  io  cm).  In  this  range,  neither  the  peak  velocity  nor  the  velocity 
spreads  1u*e  strongly  affected.  Therefore,  the  loo  tlme-of-fllght  traces 
can  be  used  as  good  indicators  of  Interactions  close  to  the  target. 

Some  energy,  albeit  a  small  amount.  Is  transferred  from  the  target 
plasma  to  the  ambient-plasma  In  the  100  mTorr  regimes  this  la  observed  In 
the  framing  camera  pictures  of  visible  light  amission.  A  fusty  shell 
expands  with  a  speed  of  about  250  km/ sec  Into  a  100  mTorr  >2  background 
gas. 


Measurements  or  the  magnetic  field  dynamics  show  that  ambient  field 
Is  largely  swept  out  of  the  region  traversed  by  the  target  plasma  and  Is 
compressed  ahead  of  It,  as  lllustrstsd  In  Fig.  10.  An  extensive' 
description  of  the  aagnetlc  rteld  behavior  is  to  be  round  In  papers  by 
Kacenjar,  et  al.3*1* 

Magnetized  ion-ion  Instability.  For  what  experimental  parameters  la 
the  magnetized  l on-lon  Instability  most  likely  to  be  observed?  To  answer 
this  question  we  examine  the  Instability  criteria  outlined  in  Lampe  et 
al.  The  two  most  stringent  instability  criteria  are  that  of  avoiding 
electromagnetic  stabilization  (Vg/V.  *H. <  2.5  )  and  fitting  at  least  one 
parallel  wavelength  in  the  system  size  <  2vR  ). 
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Figure  10.  Taaporal  behavior  of  the  change  In  aagnetic  field  detected  by 
•  Magnetic  loop  probe  located  3  ca  froa  the  target.  The 
applied  field  was  600  gauss.  Note  the  field  coapression  peak 
at  100  nsec  and  near  depletion  of  the  applied  field  (aagnetic 
bubble)  for  t  >  120  nsec.  (Aablent  gas  was  5  aTorr  hydrogen; 
laser  energy  was  T.7  J) 


Salth  and  Hubs20  quantitatively  delineated  expected  regions  of 
instability  In  laser  energy  aablent  pressure  apace,  as  illustrated  In 
Figure  11. 

TO  detect  the  presence  of  Instability,  an  array  of  four  Ion  tlae-or- 
fllght  detectors  la  deployed  perpendicular  to  the  Magnetic  field 
direction,  but  at  various  angles  froa  the  target  noraal. 

1  clear  difference  between  the  target  plaaaa  Ion  distributions  with 
and  without  the  Magnetic  rield  la  seen  In  Fig.  12.  Without  a  field 
present,  as  In  the  ion  distributions  on  the  left  aide  of  Fig.  12,  the  ion 
dlatrlbutlona  are  well-behaved  and  slallar  to  those  generated  in  a 
vacuiai.  In  contrast,  the  distributions  In  the  presence  of  the  aagnetic 
field  (on  the  right)  are  broadened  and  have  lower  velocity  (latar-tlae) 
peaks  which  are  like  the  expected  signatures  of  beaa-plasaa  instability. 
The  pressure  la  very  low  and  collisions  play  little  or  no  role.  (Me  used 
hydrogen  gas  to  have  a  low  aablent-plesaa  atoale  nuaber  In  order  to  be 
aore  susceptible  to  Instability.)  To  test  the  conjecture  that  raising  the 
ratio  of  drift  velocity  to  Alf van  velocity  above  about  2.5  quenches  the 
instability.  We  compared  Ion  distributions  resulting  froa  use  or  nitrogen 
*  4)  and  hydrogen  ■  <),  as  shown  In  Fig.  13.  The  aagnetic 

field  dependent  interaction  that  occurs  in  hydrogen  is  not  seen  in 
nitrogen.  This  lends  support  to  the  notion  that  MZI  aay  be  occurring  in 
the  carbon-hydrogen  coabl nation  but  not  In  carbon-nitrogen  under  otherwise 
Identical  conditions,  due  to  the  electromagnetic  stabilisation  criteria. 
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Figure  11.  Magnetic  ion-ion  instability  window  (qualitative),  the 

unstable  region  is  bounded  on  the  right  by  electroaagnetlc 
stabilisation  when  V-  >  2.5  V^,  on  the  left  by  the  requirement 
that  the  energy  density  in  the  target-plasma  exceed  that  of 
the  magnetic  field,  on  the  bottom  by  requiring  that  the 
transit-time  of  an  ambient  ion  passing  through  the  shell 
exceed  one  momentum  transfer  e-fold,  and  bounded  on  the  top  by 
the  practical  requirement  that  the  equal  mass  radius  fit 
inside  the  experiment.  The  top  can  also  be  limited  by  the 
cyclotron  radius  of  the  debris  ions  at  high  field  strengths. 


Finally,  sonatinas  ion  signals  are  observed  with  components  having 
higher  velocities  than  in  the  original  target-plasma  distributions,  'mis 
occurs  at  low  incident  laser  energy  in  cases  which  exhibit  the  other  signs 
of  magnetic  field  dependent  interaction.  He  have  very  few  examples  of 
these  accelerated  ions,  and,  therefore  the  results  are  considered 
tentative.  But,  these  "fast  ions"  nay  be  the  result  of  a  reflection  or 
acceleration  process  near  the  coupling  region. 

S1MURY 

The  laser-target  counterstreaaing  plasma  experiment  shows 
interactions  in  both  the  collislcnal  (high-pressure)  regime  and  the 
collisionless  (low-pressure)  regime. 

The  bulk  of  our  experimental  work  has  been  in  the  collialonal 
regime.  Here,  we  Increase  the  ambient  pressure  in  the  experiment  Into  the 
hundreds  to  thousands  of  mlllitorr  pressure  range.  Hell-formed  Taylor- 
Von-Neumann-Sedov  shocks  (or  "blast -waves'1)  are  formed  when  the  expanding 
target-plasaa  sweeps  up  the  ambient  plasma.  The  blast-wave  radius  versus 
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Figure  12.  Comparison  of  target-plasma  Ion  time- of- flight  traces  without 
a  magnetic  field  present  (left)  and  with  a  600  G  field  present 
(right).  Mote  the  distorted  debris  distributions  in  the 
presence  of  a  magnetic  field.  (Both  cases  had  hydrogen 
ambient  gas  at  15  mTorr;  shot  on  left  had  7.4  J  laser  energy 
while  that  on  right  had  fl.6  J). 


w-1 


Figure  13.  Comparison  of  target-plasma- Ion  distributions  with  atomic 
nuaber  of  the  background  gas.  (B--  600  G,  n.  •  0.01 
n0  -  6x10lz  lons/cc,  and  Vd  «  a*10'cm/ss  E.  »  8  J:  Ion 
detector  at  47*)  1 
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'.it"  dependence  follows  the  classic  Mast-wave  theory  over  a  oroao  range 
nf  operating  parameters.  The  coupling  shells  In  this  regime  also  satisfy 
other  predictions  of  blast-wave  theory  such  as  the  temperature  and  density 
jumps.  But,  the  fraction  of  the  t  -get-plasma  energy  actually  coupleo 
into  the  blast-wave  system  falls  elf  rapidly  below  about  500  sTorr 
pressure  (for  ions  with  an  Initial  speed  of  *50  km/s). 

However,  in  addition  to  smooth  classic  blast -wave  shock  regions, 
features  are  observed  which  are  not  Included  in  simple  blast-wave  models, 
ie.,  the  aneurisms,  fluff  to  the  rear  of  the  target,  and  unresolved 
turbulence  Inside  the  blast-front.  The  mechanisms  causing  this  structure 
are  still  under  active  investigation.  Other  non-blaat-wave-ltke  features 
observed  Include:  a  higher  density  Inside  the  bubble  than  expected  and  an 
electron  density  ramp  and  ledge  in  front  of  the  shook. 

another  area  in  which  very  good  progress  has  been  made  experimentally 
and  theoretically  la  in  describing  the  dynamics  of  the  magnetic  field 
under  the  influence  of  the  expanding  target-plasma  and  shock  fronts.  It 
has  been  established  that  the  uncoupled  target-plasma  compresses  some  of 
the  external  magnetic  field  ahead  of  it.  Also,  in  very  recent 
experimental  runs,  additional  compression  due  to  the  blast-wave  is 
observed. 

In  the  collisionless  regime,  the  experiments  show  magnetic  field 
dependent  interactions  between  the  fast  target  ions  and  ambient  plasma. 

In  many  ways  these  observations  are  similar  to  those  expected  from  the 
magnetized  ion-ion  instability. 

acknowledgements 

The  authors  thank  contributions  from  L.  Daniels,  V.  Salzman,  V. 
Nocerlno  and  D.  Car pen tier.  This  work  was  sponsored  by  the  Defense 
Nuclear  Agency  under  task  order  W990MXMA  and  task  order  number  00011. 


REFERENCES 

a.  Sam  Houston  State  University,  Huntsville,  TX  773*1. 

• 

1.  C.  Taylor,  Proc.  Hoy.  Soe.  A201.  159  0950)  and  1201  ■  175 

(1950). 

J.  von  Neumann  and  R.D.  Rlehmyer,  J.  Appl.  Phya.  £ ,232 

0950);  and  H.H.  Goldstlne  and  J.  von  Neumann,  Comm.  Pure  Appl. 
Math.  VIII  327  (1955). 

L.I.  Sedov.  "Similarity  and  Dimensional  Methods  in 
Mechanics,"  (Academic  Press,  NT,  1959,  ed.  M.  Holt). 

2.  B.H.  Ripln,  R.R.  Whitlock,  F.C.  Toung,  S.P.  Obenschain, 

E.A.  McLean,  and  R.  Decosta,  Phys.  Rev.  Lett.  *3,  350  (1979)i  also 
see  B.H.  Rlpin  at  ml...  Phys.  Fluids  jJi,  1012TT980)  andjj*,  990 
0981). 

3.  S.T.  Kacenjar,  B.H.  Ripln,  J.A.  Stamper,  J.  Grun,  and  E.A. 

McLean,  NRL  Memo  Report  #5260  0  98*). 

«.  J.A.  Stamper,  B.H.  Ripln,  E.A.  McLean,  and  S.P. 

Obenschain,  NRL  Memo  Report  #5278  (198*). 

5.  E.A.  McLean,  J.A.  Stamper,  H.R.  Grlem,  A.N.  All,  B.H. 

Ripln,  and  C.K.  Manka,  NRL  Memo  Report  #527*  (198*). 

6.  M.J.  Herbal  and  J.  Grun,  Phys.  Fluids  2*,  1917  (1981); 

M.J.  Her  bat  et  al.  NRL  Memo  Report  #4893  0  980)  and  #5003  (1983) 
see  also,  M.J.  Her bat  et  al.,  "Laser  Interaction  and  Related 
Plasma  Phenomena,  Vol.  6,”  pgs.  317-33*.  eds.  H.  Hora  and  C. 

Miley,  Plenua  Press  (198*). 


62 


7.  J.  Crun,  S.P  Obwnschain,  B.H.  R;;-in,  R.R.  Whitlock,  t.A. 

McLean,  J.  Gardner,  M.J.  Herb;.,  and  J.A.  Stamper,  Phys.  Fljidf  '• 
588  (1983). 

8.  J.  Grun,  R.  Decosle,  B.H.  Rlpin,  and  J.  Gardner,  Appl. 

Phys.  Leu.  j}9.  5*5  (1981). 

9.  J.  Grun,  R.  Slellingwerr,  and  B.H.  Rlptn,  to  appear  In 

Phys.  Fluids  0986). 

10.  A.W.  All  and  E.A.  McLean,  J.  Quant.  Sped.  Rad.  Trans.  33, 

381  0  985). 

D.  Y.B.  Zel'dovlch  and  l.P.  Raizer,  'Physics  of  Shock  Waves 

and  High  Temperature  Hydrodynamic  Phenomena,  Vol.  1,'  (Academic 
Press,  NT,  1966). 

12.  B.H.  Ripln,  J.A.  Stamper,  and  E.A.  McLean,  NRL  Memo  Report 

52T9  (198#). 

13.  J.A.  Stamper,  C.K.  Manka,  E.A.  McLean,  A.N.  Moslovych, 

S.P.  Obenschain,  and  B.H.  Rlpin,  (to  be  published). 

1*.  C.  Longmire,  (private  communication). 

15.  M.  Keskinen,  (private  communication). 

16.  J.  Giuliani;  NRL  Memo  Report  #5*20  (198»). 

17.  C.V  Skllzkov,  In  'Laser  Interaction  and  Related  Plasma 

Phenomena,  Vol.1"  pg.  235.  eds.  H.  Schwarz  and  H.  Kora,  Plenum 
Press,  NY  (197D. 

N.G.  Basov  et  al..  In  'Laser  Interaction  and  Related 

Plasma  Phenomena,  Vol  3B,'  pg.  553,  eds.  H.  Schwarz  and  H.  Hora. 
Plenum  Press,  NY  097#). 

D.M.  Wilke,  Los  Alamos  Report  LA-9182-T  0982). 

18.  H.  Lampe,  W.M.  Manhelmer ,  and  K.  Papadopoulos,  NRL  Memo 

Report  #3076  0  975). 

19.  S.T.  Kacenjar,  M.  Hausman,  M.  Keskinen,  A.W.  All,  J.  Grun, 

C.  Manka,  E.A.  McLean,  and  B.H.  Rlpin,  to  appear  In  Phys.  Fluids 
(1986). 

20.  R.  Smith,  and  J.  Huba,  NRL  Memo  Report  #5092  (1983). 


63 


INSTABILITY  OF  TAYLOR-SEDOV  BLAST  WAVES 
PROPAGATING  THROUGH  A  UNIFORM  GAS 


65 


Volume  66,  Number  2 1 


PHYSICAL  REVIEW  LETTERS 


27  May  1991 


Instability  of  Taylor -Sedov  Blast  Waves  Propagating  through  a  Uniform  Gas 

J.  Gran,  J.  Stamper, (t)  C.  Manka,  J.  Resnick, (b)  R.  Burris, tc*  J.  Crawford, w>  and  B.  H.  Ripin 
Space  Plasma  Branch.  Plasma  Physics  Division,  Naval  Research  Laboratory.  Washington,  D  C.  20375-5000 

(Received  10  January  1991) 

We  present  the  first  measurements  of  an  instability  in  Taylor-Sedov  blast  waves  propagating  through 
a  uniform  gas.  The  instability  occurred  in  a  gas  whose  adiabatic  index  was  low.  Amplitude  perturba¬ 
tions  grew  as  a  power  of  time.  Our  observations  are  compared  to  theory. 

PACS  numbers:  47.40.Nm.  28.70.+y.  47.20.-k,  S2.35.Tc 


Blast-wave  instability  may  contribute  to  the  structur¬ 
ing  observed  in  supernovae  and  play  a  role  in  the  forma¬ 
tion  of  stars  and  galaxies. 1  Unfortunately,  knowledge  of 
blast-wave  instabilities  is  based  almost  entirely  on 
theoretical  considerations,  and  these  have  been  accom¬ 
panied  by  considerable  controversy.2*4  Unstable  blast 
waves  had  not  been  observed  experimentally,  leading 
some  to  conclude  that  blast  waves  must  be  stable.3 

We  present  the  first  measurements  of  an  instability  in 
Taylor-Sedov  blast  waves3  propagating  through  a  uni¬ 
form  gas.  The  instability  occurred  in  an  ambient  gas 
whose  adiabatic  index  y  was  a  low  1.06  ±0.02.  Pertur¬ 
bations  grew  as  a  power  of  time.  Our  observations  are 
compared  to  a  theory  described  in  papers  by  Vishniac 
and  Ryu.4 

Blast  waves,  in  our  experiment,  are  produced  by  the 
expansion  of  ablation  plasma  from  the  surface  of  laser- 
irradiated  foils  into  an  ambient  gas  (Fig.  1).  A  6-jjm- 
thick  polystyrene  foil  is  placed  in  a  chamber  which  is 
first  evacuated  and  then  filled  to  S-torr  pressure  of  nitro¬ 
gen  or  xenon  gas.  The  foil  surface  is  heated6  to  about 
800  eV  with  a  200-J,  1.054-/im,  5 -ns  pulse  from  the 
Pharos  III  Nd-glass  laser,  which  is  focused  to  a  3- 
TW/cm2,  880-/im-diam  spot.  Ablation  plasma  from  the 
hot  foil  surface  propagates  supersonically  into  the  back¬ 
ground  gas7  at  about  700  km/s  and,  much  like  the  prod¬ 
ucts  of  a  chemical  explosion,  forms  a  blast  wave. 
(Simultaneously,  the  background  gas  is  photoionized  by 
radiation  from  the  vicinity  of  the  laser’s  focal  point.) 


FIG.  1.  Experimental  setup. 


We  have  verified  through  extensive  experimentation  that 
this  laser-ablation  method  forms  classical  Taylor-Sedov 
blast  waves  when  the  interaction  between  the  ablation 
plasma  and  the  gas  is  collisional,  and  when  the  mass  of 
the  swept-up  ambient  gas  is  greater  than  the  mass  of  the 
ablation  plasma.*  For  nitrogen,  collisional  coupling 
occurs  at  gas  pressures  exceeding  0.3  torr. 

Blast-front  structure  is  photographed  using  the  well- 
known  dark-field  imaging  method,  which  is  sensitive  to 
the  square  of  fluctuations  in  the  index  of  refraction.  In 
our  implementation,  a  0.53-/im,  <  1  -ns-duration,  5-cm- 
diam,  few-mJ  laser  probe  illuminates  the  blast  wave 
side-on.  Electron-density  gradients  within  the  blast  wave 
deflect  a  part  of  the  probe  while  the  remainder  passes 
through  undisturbed.  The  probe  beam  emerging  from 
the  blast-wave  region  is  then  relayed  onto  a  film  surface 
with  a  telescope.  A  stop  placed  at  a  focal  point  inside 
the  telescope  blocks  the  undisturbed  component  of  the 
probe  light  but  passes  the  deflected  part,  thereby  forming 
an  image  in  wh::h  fluctuations  in  the  index  of  refraction 
(and  hence  electron  density)  appear  as  bright  features  on 
a  dark  background. 

In  addition,  visible  emission  from  the  blast  front  is 
photographed  with  a  very  fast  (120-ps  to  5-ns  gate 
time),  four-frame,  microchannel-plate  intensifier  camera. 
This,  together  with  the  dark-field  image,  provides  five 
photographs  per  shot  of  the  blast  wave  at  different  times 
in  its  evolution.  The  dynamics  of  the  blast  wave  are 
reconstructed  by  combining  the  results  of  dark-field  and 
emission  photographs  taken  at  different  times  on  indivi¬ 
dual  and  multiple  shots.  Also,  the  spectrum  of  light  em¬ 
itted  by  the  ambient  gas  before  and  after  the  passage  of 
the  blast  wave  is  measured  with  temporal  and  spatial 
resolution.  To  do  this,  we  image  a  3-mm-diam  spot  in 
front  of  the  foil  onto  the  slit  of  a  spectrometer  and  record 
the  resulting  spectrum  with  a  streak  camera. 

We  find  that  Taylor-Sedov  blast  waves  formed  in  ni¬ 
trogen  gas  are  always  stable  and  smooth — like  the  exam¬ 
ple  in  Fig.  2(a).  In  startling  contrast  [Fig.  2(b)],  the 
surface  of  the  blast  fronts  launched  into  an  ambient  xe¬ 
non  gas  is  wrinkled  like  a  dried  prune.  This  wrinkling  is 
quantified  by  tracing  and  then  Fourier  transforming9  the 
outer  edge  of  the  front,  which  is  equivalent  to  looking  at 
the  projection  onto  a  plane  of  the  edge  of  an  unstable 
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FIG.  2.  Dark-field  shadowgraph  of  (a)  a  stable  blast  wave  in  nitrogen  gas,  and  (b)  an  unstable  blast  wave  in  xenon  gas  (at 
243  ns). 


sphere.  The  results  are  presented  as  A,(k)/R  vs 
logio(£R),  where  A,(k)  is  the  full  amplitude  of  the 
mode  with  wave  number  k  at  time  r,  and  R  is  the  aver¬ 
age  radius  of  the  blast-front  boundary. 

The  blast-wave  trajectory  in  xenon  and  the  evolution 
of  instability  amplitude  A,ik)/R  are  shown  in  Fig.  3. 
From  6  to  18  ns  the  front  moves  at  a  constant  speed  cor¬ 
responding  to  the  velocity  of  the  ablation  r’asma,  a  blast 
wave  not  yet  having  formed.  The  front  is  slightly  struc¬ 
tured,  but  these  nonuniformities  do  not  grow.  By  23  ns  a 
blast  wave  propagating  with  the  tVi  Taylor-Sedov 
dependence  has  formed.  Now  the  surface  becomes 
significantly  more  wrinkled  and  spikelike  protuberances 
shoot  ahead  of  the  front:  A,{k)(R  increases  as  a  power 
of  time  until  300  ns.  It  is  noted  that  the  blast  wave  does 
not  fall  apart  or  otherwise  dissipate,  but  propagates  as  a 
shocklike,  albeit  structured,  front.  As  the  protuberances 
get  larger  they  become  increasingly  more  difficult  to  ob¬ 
serve.  By  400  ns  they  are  not  seen  at  all  and  the  blast 
wave  takes  on  the  appearance  of  a  slightly  structured  but 
basically  stable  shock. 

A  power  law  of  the  form  •4/(fc)//?5K/s(**)  was  fitted 
to  the  ^,(&)/R-vs-time  data  during  the  period  of  growth, 
with  the  results  shown  in  Fig.  3(c).  We  find  clear 
growth  for  modes  satisfying  0.7  <  In (kR)  <  2.  Max¬ 
imum  growth  occurs  at  ln(k/?)"l,  where  S “1.6,  and 
minimum  growth,  with  S”  0.3,  occurs  at  In  (kR)  •2. 
The  fit  by  a  power  law  is  very  good  (correlation  >0.7) 
for  In(kR)  <  l.S,  but  worse  (i.e.,  data  are  more  noisy) 
for  larger  values  of  In  (kR).  Noise  may  be  the  reason 
why  S(kR )  stays  clamped  at  0.3  for  2  <  In  (kR)  <  3  and 
does  not  decrease  to  zero. 

A  basic  difference  between  stable  blast  waves  in  nitro¬ 
gen  and  unstable  blast  waves  in  xenon  is  that  the  former 
propagate  in  a  gas  with  adiabatic  index  7n"“  1.3  ±0.1, 
and  the  latter  in  a  gas  with  yXe  “  1 .06  ±  0.02.  To  derive 
Txe  we  utilize  the  observation  that  blast  waves  in  nitro¬ 


gen  and  xenon  propagate  according  to  the  Taylor-Sedov 
blast-wave  relation 10 


dec 


7S(y— l)(y+l)2 
16jr(3y— l)po 
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,2/5 


where  d  is  the  distance  between  the  focal  spot  and  the 
blast-wave  front  and  po  is  the  gas  mass  density.  There¬ 
fore,  by  dividing  the  measured  d  in  nitrogen  by  the  mea¬ 
sured  d  in  xenon  at  any  given  time  we  arrive  at  a  rela¬ 
tionship  between  yXe,  7n,  and  the  mass  of  each  gas 
species.  Solving  numerically  for  yXc  as  a  function  of  7n, 
we  find  that  as  7n  varies  from  1  to  f ,  7Xc  varies  from  I 
to  1.13.  Hence,  for  any  reasonable  value  of  7n  the  value 
of  yXc  is  less  than  1.13.  In  past  experiments8  we  have 
measured  7n  to  be  1.3  ±0.1,  which  implies  that 

7Xe- 106  ±0.02. 

7Xc  is  lower  than  7n  because  prior  to  the  arrival  of  the 
blast  wave  xenon  gas  radiates  much  more  than  nitrogen 
gas.  (Radiation  increases  the  degrees  of  freedom  within 
a  gas  and  hence  reduces  its  effective  7. )  This  is  demon¬ 
strated  by  examining  the  spectrum  of  light  emanating 
from  a  spot  in  front  of  the  laser's  focal  point  [Fig.  3(d)]. 
In  nitrogen  gas  the  laser-induced  explosion  produces  lit¬ 
tle  measurable  emission  prior  to  the  arrival  of  the  blast 
wave.  Immediately  before  the  blast  wave  arrives  at  the 
observation  point  there  is  a  slight  increase  of  Nl+  and 
N2+  lines:  These  lines  are  probably  excited  by  UV  or 
heat  from  the  blast  wave.  When  the  blast  front  arrives, 
there  is  a  sudden  increase  in  emission  from  the  N l+, 
N2+,  and  target  C2+  lines,  as  well  as  an  increase  in  con¬ 
tinuum  emission.  In  contrast,  xenon  emits  copiously  in 
many  Xel+  and  some  Xe2+  lines  from  the  moment  the 
laser  strikes  the  foil.  Arrival  of  the  unstable  blast  wave 
is  signaled  by  a  more  gradual  increase  in  continuum 
emission,  but  line  emission  is  not  changed  significantly. 
We  conclude,  therefore,  that  it  is  the  radiation  in  xenon 
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FIG.  3.  Instability  growth,  (a)  Blast-wave  trajectory,  (b)  Amplitude  growth  as  a  function  of  time  for  different  values  of  In(fcA). 
The  scatter  of  points  between  30  and  43  ns  at  each  value  of  ln(£A)  indicates  the  shot-to-shot  reproducibility  in  this  experiment,  (c) 
Growth  exponent  S(kR).  (d)  Emission  spectra  as  a  function  of  time  in  nitrogen  and  xenon  gas  from  a  spot  which  is  2  cm  from  the 
laser’s  focal  point. 


which  reduces  its  effective  y  below  that  of  nitrogen. 

Among  the  various  theories,2-4  the  one  by  Vishniac 
and  Ryu4  predicts  that  Taylor-Sedov  blast  waves  in  a 
uniform  gas  with  y  <  1 .2  will  become  unstable.  In 
spherical  geometry,  perturbations  are  predicted  to  grow 
as  Yim(6,p)ts{,\  where  Y/„  are  spherical  harmonic 
modes.  In  planar  geometry,  perturbations  grow  as 
etky{sikx)'  wf,cre  x  is  the  direction  of  blast  propagation. 
For  y— 1.1,  numerically  calculated  maximum  growth 
occurs  at  ln(/)  ss  1 .5  with  RetS(/)l  as0.5,  and  ln(Lx) 
a  1.2  with  RelS(Lx)]a0.3,  for  the  spherical  and  pla¬ 
nar  cases,  respectively.  For  y*1.06,  a  less  precise  but 
analytic  calculation"  in  spherical  geometry  gives  max¬ 
imum  growth  at  ln(/)sl.7,  with  Re{S(/)]s0.7.  Theo¬ 
ry  does  not  treat  the  large-amplitude  regime  where  satu¬ 
ration  or  stabilization  may  occur. 

The  basic  predictions  of  this  theory  agree  with  our  ob¬ 
servations  of  ts  growth  in  a  low-y  uniform  gas,  but  there 
are  differences  in  the  details.  For  y— 1.06  we  measure 


maximum  growth  at  ln(fc/?)ssl  with  S&1.6,  whereas 
theory  predicts  maximum  growth  at  higher  ln(/)  and 
lower  S:  Theory  predicts  virtually  no  growth  (S  as  0.04) 
at  ln(/)  “  1.  Possible  stabilization  is  seen  at  400  ns.  We 
remind  the  reader,  however,  that  the  wave  numbers  in 
our  experiment  were  obtained  from  the  projection  onto  a 
plane  of  the  edge  of  an  unstable  sphere.  Thus,  kR  is  not 
identical  to  either  /  or  kx.  Also,  there  are  a  few  alterna¬ 
tive  explanations  for  the  observed  apparent  stabilization. 
One  may  argue,  for  example,  that  the  stabilization  is  a 
real  nonlinear  phenomenon;  or  that  y  is  large  at  some 
point  far  from  the  focal  spot  and  the  blast  wave  stabi¬ 
lized  because  it  moved  into  that  region;  or  that  the  insta¬ 
bility  is  still  growing  but  the  density  gradients  in  the 
spikes  are  too  gentle  to  measure  at  late  times.  These  is¬ 
sues  will  be  addressed  in  future  work. 

We  point  out  that  the  phenomenon  described  here  is 
not  the  Rayleigh-Taylor  instability  commonly  associated 
with  decelerating  systems.  That  instability  is  caused  by 
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opposing  density  and  pressure  gradients  and  exhibits  ex¬ 
ponential  growth.  This  instability  is  associated  with  the 
sloshing  of  material  within  the  blast-wave  shell  and  ex¬ 
hibits  power-law  growth.  Our  instability  mechanism 
may  be  briefly  explained  as  follows:4  The  thermal  pres¬ 
sure  which  drives  the  blast  wave  is  perpendicular  to  the 
local  blast  shell  surface,  while  the  external  ram  pressure 
(ambient  density  times  square  of  blast  speed)  is  antipar¬ 
allel  to  the  direction  of  propagation.  In  a  uniform  blast 
wave  the  front  and  the  direction  of  propagation  are  or¬ 
thogonal  so  that  the  thermal  and  ram  pressures  are  both 
perpendicular  to  the  shell  surface.  But  in  a  rippled  blast 
wave  the  thermal  pressure  is  no  longer  parallel  to  the 
propagation  direction  while  the  ram-pressure  orientation 
does  not  change.  Therefore,  there  appears  a  net  pressure 
along  the  blast  front  surface  that  accelerates  mass  into 
the  lagging  trough  parts  of  the  ripple.  The  now  heavier 
trough  has  more  momentum,  slows  less,  and  consequent¬ 
ly  moves  ahead.  Then  the  process  reverses  and  an  oscil¬ 
lation  ensues.  This  oscillation  is  damped  out  at  wave¬ 
lengths  comparable  to  the  blast  shell  radius  and  wave¬ 
lengths  much  larger  than  the  shell  thickness.  Otherwise, 
it  grows.  Since  shell  thickness  is  related  to  y,  so  is  the 
growth  exponent  of  the  instability. 

In  conclusion,  we  have  shown  experimentally,  for  the 
first  time,  that  Taylor-Sedov  blast  waves  in  a  uniform 
gas  are  unstable  when  the  adiabatic  index  y  is  suf¬ 
ficiently  low,  in  our  case  1.06.  Perturbed  amplitudes 
grow  as  a  power  of  time.  Our  results  confirm  the  basic 
predictions  of  Vishniac  and  Ryu. 
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and  Dr.  David  Book  for  reviewing  early  versions  of  this 
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Observation  of  high-pressure  blast-wave  decursors 
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We  observe,  for  the  first  time  in  a  laboratory,  the  formation  of  decursors  at  intersections  of 
planar  surfaces  with  powerful  blast  waves.  The  blast  waves,  which  have 
hundreds-kilobar  overpressures,  are  created  by  rapid  ablation  of  material  heated  with  an 
intense  laser  beam. 


It  is  well  known  that  a  blast  wave  colliding  with  a 
planar  surface  sets  up  a  reflection  that  propagates  back  into 
the  ambient  gas  heated  by  the  passage  of  the  incident  front. 
The  interacting  incident  and  reflected  blast  fronts  can  form 
structures  known  as  Mach  stems  and  triple  points.1  Even 
more  complex  structuring  can  occur  if  the  reflecting  sur¬ 
face  is  shaped,  or  if  the  sound  speed  near  the  surface  varies. 
For  example,  if  the  sound  speed  near  the  reflecting  surface 
is  higher  than  in  the  rest  of  the  ambient  gas,  part  of  the 
blast  front  moves  faster  and  outruns  the  bulk  of  the  blast 
wave.  This,  in  turn,  forms  additional  shocks  and  compli¬ 
cated  vortical  flow  patterns.  The  structure  which  outruns 
the  blast  wave,  called  a  precursor,  has  been  studied  at 
length  both  theoretically  and  experimentally.2’3  It  is  also 
possible  to  envision  situations  in  which  the  sound  speed 
near  the  reflecting  surface  is  slower  than  in  the  rest  of  the 
ambient  gas.  Under  such  circumstances,  one  expects  that 
one  part  of  the  blast  wave  will  lag  behind  the  rest.  The 
existence  of  such  structure,  called  a  decursor,  has  been 
considered  theoretically  but  it  has  never  been  seen  in  an 
experiment. 

We  report  the  first  observation  of  blast-wave  decursors. 
In  our  experiment  Mast  waves  are  created  by  rapid  ablation 
of  material  heated  with  an  intense  laser  beam.  We  will 
show  that  such  laser-induced  blast  waves  follow  Taylor- 
von  Neumann-Sedov  scaling4-3  and  that  their  overpres¬ 
sures  are  much  higher  than  overpressures  achievable  with 
chemical  explosives. 

Figure  1  shows  how  a  laser  is  used  to  create  blast 
waves.  A  small,  millimeter-diameter,  foil  or  sphere  target  is 
placed  in  an  evacuated  chamber,  which  is  then  backfilled 
with  gas.  The  target  surface  is  irradiated  at  high  intensity 
by  a  focused  pulse  from  the  PHAROS  III  Nd  glass  laser 
(which  has  three  beams  and  can  produce  up  to  1500  J  per 
pulse).  Material  ablating  from  the  target  surface  expands 
rapidly  into  the  background  gas — much  like  the  products 
of  a  chemical  explosion — and  forms  a  blast  wave.4,5  In  the 
experiments  described  here,  the  targets  are  9-/zm-thick 
polystyrene  foils  and  1/16  in.  diam  aluminum  rods,  the 
atmosphere  is  5  Torr  of  nitrogen  gas,  and  the  laser  pulse  is 
a  100-200  J,  1.054  pm,  5  ns,  pulse  focused  to  a  880-jtm- 
diatn  spot  This  produces  an  irradiance  of  3  TW/cm2.  With 
these  parameters  the  foil  surface  heats  to  about  800  eV,* 
and  ablates  (or  explodes)  away  at  about  700  lcm/s.7'9  A 
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classic  blast  wave  forms  within  about  10  ns. 

The  blast  front  structure  is  photographed  using  the 
dark-field  imaging  method,  which  is  sensitive  to  density 
gradient  fluctuations.10  In  our  implementation,  a  0.53  pm, 
<  1  ns  duration,  2-in.-diam,  few-millijouk  laser  probe  illu¬ 
minates  the  blast-wave  side-on.  Density  gradients  within 
the  blast  wave  deflect  part  of  the  probe  while  the  remainder 
passes  through  undisturbed.  The  emerging  probe  beam  is 
then  relayed  onto  a  film  surface  with  a  telescope.  A  stop 
placed  at  a  focal  point  inside  the  telescope  blocks  the  un¬ 
disturbed  component  of  the  probe  light  but  passes  the  de¬ 
flected  part,  thereby  forming  an  image  in  which  blast-front 
gradients  appear  as  bright  features  on  a  dark  background. 

In  addition,  the  visible  emission  from  the  blast  front  is 
photographed  with  a  very  fast  (100  ps  shutter),  four- 
frame,  microchannel-plate  intensifier  camera.  Thus,  to¬ 
gether  with  the  dark-field  image,  we  get  five  photographs  of 
the  blast  wave  at  different  times  in  its  evolution.  This  min¬ 
imizes  shot-to-shot  variations  and  gives  an  accurate  mea¬ 
sure  of  the  blast-wave  region  history. 

Different  blast-wave-surface  configurations  are  shown 
in  Fig.  2.  Fipure  2(a)  shows  a  blast  wave  soon  (185  ns) 
after  it  hit  an  aluminum  plane  placed  6  mm  below  the 
laser’s  point-of-focus.  The  blast  wave  shows  up  as  a  bright 
circle  surrounding  the  focus.  Above  the  aluminum  plane 
are  seen  diffuse,  bright  features  that  rise  to  the  point-of- 
focus  immediately  below  the  explosion,  but  stay  closer  to 
the  plane  as  the  lateral  distance  from  the  explosion  in- 
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FIG.  1.  Experimental  setup.  A  blast  front  is  launched  by  an  energetic, 
laser-produced  plasma  expanding  into  a  low-pressure  ambient  gas.  Emis¬ 
sion  from  the  blast  front  is  imaged  with  a  four-frame,  fast-gated  intensi¬ 
fier,  optical  imager  (GOI).  Density  gradients  are  recorded  with  a  dark- 
field  shadowgraphy  camera. 
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FIG.  2.  Structuring  of  bint  warn  near  reflecting,  aluminum  surfaces  as 
seen  witb  dark-field  shadowgraphs  and  visible  emission  photographs.  The 
Mast  waves  originate  about  6  mm  above  the  surface  and  propagate  is  S 
Torr  of  nitrogen,  (a)  Spherical  blast  wave  originating  at  the  tip  of  an 
aluminum  rod  as  observed  185  ns  after  explosion-  The  laser  energy  was 
210  I.  (b)  Mach  stem,  triple  point,  and  decursor  formation  is  a  blast 
wave  originating  at  the  tip  of  an  aluminum  rod.  The  laser  energy  was  100 
J  and  the  observation  was  at  302  ns.  (c)  A  blast  wave  from  a  plastic  target 
shows  Mach  stems,  triple  points,  and  strong  decursor  at  357  os  after 
explosion.  Laser  energy  was  103  J.  (d)  Goseups  of  different  decursor 
shapes  observed  in  this  experiment 

creases.  Similar  features,  which  surround  the  point-of-fo- 
cus  itself  are  seen  in  both  the  dark-field  and  emission  pho¬ 
tographs.  These  diffuse  features,  which  do  not  exhibit  any 
sharp  characteristics  associated  with  blast  waves  and 
shocks,  are  signatures  of  material  ejected  from  the  target 
and  also  of  material  ablated  or  vaporized  from  the  surface 
of  the  reflecting  plane  by  radiation  from  the  focus  region  or 
the  blast  front  The  electron  density  of  the  ionized  material 
near  the  plane  is  comparable  to  the  density  of  the  S  Torr 
background  nitrogen  gas — implying  that  the  gas  above  the 
plane  is  mostly  singly  ionized.  (Electron  density  was  esti¬ 
mated  by  calculating  the  minimum  probe  deflection  neces¬ 
sary  to  miss  the  stop  inside  the  dark-field  telescope. )  The 
reflected  part  of  the  blast  wave  is  buried  within  the  features 
of  the  ionized  gas  and  is,  therefore,  not  clearly  visible.  A 
dark  circle  near  the  top  of  the  photograph  marks  the  pas¬ 
sage  of  the  laser  beam. 

Figure  2(b)  taken  later  at  302  ns,  shows  well  devel¬ 
oped  Mach  stems  and  triple  points.  As  in  the  previous 
picture,  there  is  evidence  of  ionized  material  around  the 
focal  point  and  near  the  reflecting  plane.  Where  the  blast 
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wave  meets  the  plane  we  see  the  rounded,  tucked-back 
features  of  a  decursor. 

Much  clearer  and  stronger  decursors  are  seen  in  Fig. 
2(c),  observed  at  357  ns  after  explosion.  By  this  time  the 
diffuse  plasma  features  near  the  plane  are  no  longer  seen. 
The  decursor  is  also  straighter  than  the  rounded  decursor 
in  the  previous  shot.  The  same  photograph  also  shows 
Mach  stems  and  triple  points. 

A  likely  explanation  of  why  our  experiment  forms  de¬ 
cursors  is  that  the  material  ablated  or  vaporized  from  the 
reflecting  plane  forms  a  layer  with  increased  density 
and/or  effective  mass,  and  a  lower  sound  speed.  The  source 
of  energy  for  this  ablation  is  radiation  from  the  laser- 
heated  plasma  and/or  the  blast  wave.  Radiant  energy  from 
plastic  and  aluminum  targets  which  are  heated  to  800  eV 
constitute  approximately  5%  and  25%  of  the  incident  laser 
energy,  respectively.11  The  ratio  of  blast-wave  speed  to  de- 
cursor  speed  can  be  estimated  by  assuming  that  it  is 
roughly  the  same  as  the  ratio  of  their  respective  distances 
from  the  point  of  explosion.  For  the  shots  shown  in  Fig.  2 
that  ratio  is  1.1. 

The  spatial  and  temporal  scales  of  laser-induced  blast 
waves  are  much  shorter  than  the  scales  in  traditional  shock 
tube  or  explosively  driven  experiments.  Typical  laser-labo¬ 
ratory  scales  are  centimeter  and  submicrosecond  versus 
tens  of  centimeters  and  milliseconds  in  the  more  traditional 
methods.  Also,  our  experiments  were  done  at  a  fraction  of 
atmospheric  pressure.  It  is,  therefore,  legitimate  to  ask  how 
laser  driven  shocks  relate  to  those  produced  in  more  usual 
environments.  We  proceed  to  show  that  laser-produced 
blast  waves  follow  the  scaling  laws  of  ideal  blast  waves.  We 
will  also  show  that  laser-produced  blast  waves  achieve 
overpressures  unreachable  by  explosively  driven  methods. 

Ideal  blast  waves  propagate  according  to  the  Taylor- 
von  Neumann-Sedov  relation:12 

R*f(ff/p0),/5/2/s, 

with 

f={[75 (r-  l)(r+  l)2]/16ir(3r-  1)},/5.  (1) 

where  E  is  the  explosion  energy,  p$  is  the  ambient-gas  den¬ 
sity,  t  is  the  time  of  observation,  and  y  is  the  ratio  of 
specific  heats.  These  expressions  were  derived  assuming 
spherical  symmetry,  energy  and  momentum  conservation, 
and  an  instantaneous,  massless,  point  explosion.  We  have 
performed  detailed  experiments  to  verify  that  laser-pro¬ 
duced  blast  waves  follow  Eq.  ( 1  ).4  Many  parameters  such 
as  the  laser  energy,  laser  focal-spot  size  (and  therefore  the 
ablation  velocity),  pressure,  and  observation  times  were 
varied  in  these  experiments.  As  long  as  the  ambient  pres¬ 
sure  is  high  enough  (  >  0.5  Torr)  to  place  the  experiment 
in  a  collisional  regime,  laser-produced  blast  waves  follow 
the  Taylor-von  Neumann-Sedov  relation.  From  the  data 
and  Eq.  (1)  we  determined  that  f  =  1.0  ±0.1  in  a  nitrogen 
ambient  gas,  which  implies  a  y  of  1.4  ±0.2. 

Blast-wave  overpressure  is  determined  from  the  dark- 
field  photographic  data  using  the  following  expressions:212 
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FIO.  3.  Later-created  blast-wave  overpressures  compered  to  ideal  over- 
premnres  tad  overpressures  created  by  chemical  explosives.  Blast  waves 
created  by  nuclear  explosions  follow  the  ideal  curve.  The  horizontal  axis 
it  the  distance  from  the  explosion  point  scaled  with  (ambient  pressure/ 
yield)171. 
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where  Pis  the  peek  overpressure,  P0  is  the  ambient  pres¬ 
sure,  M  is  the  blast-wave’s  Mach  number,  and  Q  is  the 
sound  speed  in  the  ambient  gas. 

Figure  3  plots  the  overpressures  inferred  from  Eqs. 
(2)-(4)  as  a  function  of  distance  from  the  explosion  scaled 
using  the  Sachs  scaling  law.1,13  This  law,  which  is  based  on 
geometrical  similarity  principles  and  conservation  of  mo¬ 
mentum,  is  used  to  plot  and  compare  experiments  per¬ 
formed  with  different  explosive  yields  and  at  different  am¬ 
bient  pressures.  Also  plotted  are  curves  expected  of  an 
ideal  Mast  wave  and  a  blast  wave  created  with  explosives.’ 
This  figure  shows  that  laser-produced  blast  waves  follow 
the  ideal  scaling  law,  as  do  nuclear  blast  waves2  which 
occur  at  higher  ambient  pressures  and  have  significantly 
larger  spatial  and  temporal  scales.  They  also  continue  to 
follow  the  ideal  law  at  much  higher  pressures  than  blast 
waves  produced  with  chemical  explosives.  Furthermore, 


the  laser  experiment  produces  blast  waves  with  overpres¬ 
sures  greater  than  100  Irbar,  which  is  more  than  1000  times 
greater  than  the  capability  of  conventional  explosives. 

In  conclusion,  we  have  shown  that  rapid  ablation  pro¬ 
duced  with  focused  laser  beams  can  serve  as  a  useful  driver 
for  the  study  of  high-pressure  blast  waves.  In  particular, 
we  demonstrated  that  such  “explosions”  produce  Taylor- 
von  Neumann-Sedov  blast  fronts  in  a  pressure  regime 
much  higher  than  otherwise  achievable  in  a  laboratory. 
These  laser-driven  blast  waves  can  be  used  to  simulate  and 
investigate  many  interesting  and  relevant  blast  front-sur¬ 
face  interactions.  We  report,  for  the  first  time,  observations 
of  blast-wave  decursor  formation.  In  future  experiments  we 
will  study  the  phenomena  observed  here  more  quantita¬ 
tively  and  will  extend  our  measurements  to  higher  ambient 
pressures. 
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INTRODUCTION 


Most  flows  in  nature  are  turbulent.  Often,  turbulence  critically  influences  die*- mire 
of  an  important  event.  For  example,  satellite  communications  can  be  interfered  with  oy 
natural  or  man-made  turbulence  in  the  earth's  atmosphere';  x  and  gamma  rays  in  the  1987A 
supernova  appeared  prematurely  because  of  turbulence1;  the  rate  of  star  formation  depends 
on  the  nature  of  interstellar  turbulence1;  turbulent  mixing  influences  the  efficiency  of 
inertial-confinement-fusion  pellets*.  There  are  many  other  examples  in  aeronautics, 
chemistry,  and  combustion. 


Turbulence  has  been  studied  for  more  than  a  century  but  h  is  still  incompletely 
understood1.  This  is  not  surprising:  Turbulence  is  by  definition  highly  nonlinear, 
experiments  are  difficult  to  interpret  quantitatively,  and  it  is  hard  to  relate  experimental 
observations  to  theoretical  predictions.  Many  experiments,  such  as  those  in  shock  tubes, 
suffer  from  extraneous  effects  caused  by  walls  and  membranes.  However,  with  rapid 
advances  in  computer  architecture,  advances  in  computer  aided  visualization  and  analysis, 
and  advances  in  experimental  methods,  rapid  progress  in  understanding  turbulence  is  now 
possible. 


Ltter  Inunction  md  Reiaied  Mount  Phenomena,  Voi.  9 
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High-power  lasers,  such  as  those  used  for  inertial -confinement-fusion  research,  are 
a  new  and  very  good  tool  for  studying  turbulence  and  other  hydrodynamic  phenomena*14 
Laser  hydrodynamics  experiments  are  not  hampered  by  walls  and  membranes  that  make 
clean  shock  tube  experiments  difficult.  Areas  of  parameter  space  not  accessible  by 
conventional  methods  can  be  reached  easily  with  high-power  lasers.  For  example, 
multi-megabar  pressures,  Mach  numbers  of  a  few  hundred,  temperatures  of  many  electron 
volts  can  be  achieved.  Intense  x-ray  radiation  can  be  turned  on  when  needed.  Also,  it  is 
relatively  easy  to  vary  parameters,  such  as  pressure,  flow  speed,  Mach  number, 
temperature,  fluid  composition,  etc.,  over  a  broad  range.  Laser  experiments  can  duplicate 
many  astrophysical  parameters* .  Also,  many  of  the  sophisticated  diagnostics  which  have 
been  developed  over  the  past  twenty  five  years  for  fusion  research  can  be  utilized  in 
hydrodynamics  experiments. 


We  use  the  N.R.L.  Pharos  HI  laser  and  target  facility  to  study  linear  and  nonlinear 
hydrodynamic  flows.  This  paper  describes  our  experimental  work  in  very  high  Mach 
number  turbulence. 


DESCRIPTION  OF  EXPERIMENT 


The  Pharos  m  laser  facility  10  has  a  three-beam  glass  laser  with  an  output  of  500 
joules  per  beam  at  a  wavelength  of  1.054pm  in  a  2-6ns  pulse.  The  laser  wavelength  can 
also  be  halved.  Imdiances  of  300  terawatts  per  square  centimeter  can  be  produced  by 
focusing  the  laser  a  few  hundred  micron  diameter  spot.  The  spatial  irradiance  profile  can  be 
stuped  and  smoothed  to  a  few  percent  level  with  the  induced-spatial-incoherence  (LSI) 
method".  Experiments  are  performed  in  a  large,  80  cm  radius  chamber  which  can  be 
evacuated  or  filled  with  gas.  A  magnetic  field  of  up  to  10  kgauss  over  more  than  100  cm* 
surrounding  the  point  of  laser  focus  can  be  turned  on.  Figure  1  shows  a  photograph  of  the 
area  around  the  Pharos  lU  target  chamber. 


FIGURE  1 .  Photograph  of  the  area  around  the  Pharos  HI  target  chamber.  Laser 

focusing  optics  and  parts  of  laser  beam  diagnostics  are  seen  in  the 
background. 
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In  litis  work  one  beam  of  the  Inset  is  focused  onto  the  surface  of  thin.  3-2Qpm. 
polvstyrcne  (CH )  or  aluminium  foils  to  produce  an  irradtance  of  1-300  terawatts/cm:  m  a 
Kid-  lOOOptn  diameter  spot  The  chamber  is  filled  with  pas  but  an  external  magnetic  field 
is  not  used.  This  is  what  happens  when  the  laser  irradiates  the  foil:11  The  irradiated  pan  of 
the  foil  (target)  is  heated  on  its  surface  to  a  temperature  of  about  800  eV'\  causing  the 
surface  to  ablate.  The  ablated  material  expands  toward  the  laser  at  about  700  kilometers 
per  second  in  a  narrow  shell14.  Half  of  the  ablation  mass  is  contained  within  a  40  degree 
half-cone  angle”  ”.  Meanwhile,  the  background  gas  is  photoionized  by  radiation  from  the 
target  or  by  impact  with  ablation  material.  A  Taylor- Sedov  blastwave  forms  in  front  of  the 
target  when  the  ablated  material  sweeps  up  a  gas  mass  equal  to  a  few  times  its  own”.  The 
non-ablated  pan  of  the  target  is  accelerated  away  from  the  laser  by  the  rocket-like  effect  of 
the  ablation  material  to  speeds  of  10-100  km/sec.  During  acceleration  the  target  heats  to 
1-2  eV  and  becomes  Rayleigh-Taylor  unstable *1*.  Later,  the  unstable  target  becomes 
turbulent  and  entrains  background  gas  so  that  it  too  becomes  turbulent.  When  the 
acceleration  ceases  the  target  decompresses  thermally.  An  example  of  a  target  irradiated  in 
a  background  gas,  showing  both  the  blast  wave  and  the  turbulence,  is  shown  in  Fig.2. 


FIGURE  2.  A  dual-time  Schlieren  shadowgraph  of  an  aluminium  foil  irradiated 
with  36  joules  in  3  ton  of  nitrogen.  The  two  exposures  were  at  33ns 
and  160ns.  Note  the  blastwave  on  the  laser  side  of  the  foil  and  the 
turbulent  region  behind  the  target. 


Our  most  important  diagnostic  tool  is  a  laser  probe  (0.53 pm,  350psec  FWHM) 
which  is  transmitted  through  die  turbulent  region  as  shown  in  figure  3.  The  phase  and 
amplitude  of  the  probe  may  both,  in  theory,  be  affected  by  the  turbulence.  The  disturbed 
probe  light  is  recorded  using  various  shadow graphy  methods,  such  as  bright-field, 
Schlieren,  or  phase-contrast  photography,  and  with  holographic  interferometry.  In 
addition,  visible  light  emission  from  the  turbulent  region  is  photographed  with  a  gated 
optical  imager  (shuttered  at  600psec)  which  is  filtered  so  that  h  images  emission  from  only 
the  target  material,  only  the  background  gas,  or  both.  Emission  spectra  measured  with  a 
temporally  and  spatially  resolved  optical  multichannel  analyzer  identify  ionic  species  and 
help  us  estimate  temperature. 
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FIGURE  3.  Schematic  of  experiment  and  main  diagnostics. 


CHARACTERISTICS  OF  THE  TURBULENT  REGION 


In  this  section  we  address  the  following  questions  about  the  turbulent  region:  How 
does  it  move  7;  What  is  turbulent,  the  target,  the  background  gas  or  both  7;  What  is  its 
spatial  wavelength  spectrum?;  and.  Can  the  spectrum  be  manipulated? 

How  does  It  move?  Figure  4  shows  a  time  sequence  of  phase-contrast 

photographs  of  the  turbulent  region.  The  speed  of  the  turbulent  front  is  determined  from 
plots  of  the  distance  between  the  original  foil  position  and  the  tip  of  the  turbulent  front 
versus  the  time  of  the  observation  (Fig.  3).  From  such  time  sequences  we  see  that  the  front 
moves  at  a  speed  which  is  somewhat  larget(<  20%)  than  thespeed  we  calculate  for  the 
target  material  itself.  Furthermore,  the  front  is  not  slowed  down  by  collisional  interaction 
with  the  background  gas.  The  front  speed  in  figure  5  is  moving  at  about  Mach  70  with 
respect  to  the  unheated  background  gas.  Mach  numbers  of  up  to  200  were  observed  in 
other  cases.  The  speed  of  the  turbulent  region  does  not  depend  strongly  on  the  background 
gas  pressure. 

What  is  turbulent?  The  fact  that  turbulence  exists  -  and  yet  its  motion  is 

not  perceptively  slowed  down  by  the  background  gas  •  may  indicate  that  a  gas-target 
interchange  instability  is  not  a  dominant  factor  in  its  development  Rather,  the  target  itself, 
which  we  know  to  be  linearly  unstable*'74,  may  be  the  primary  cause  of  the  turbulence. 

This  would  explain  why  the  turbulent  front  speed  is  somewhat  higher  than  the  target  speed 
and  why  it  does  not  depend  significantly  on  the  background  gas  pressure.  The  front  is 
probably  a  shock  driven  slightly  ahead  of  the  target  material  •  much  like  a  shock  moving 
ahead  of  a  supersonic  piston.  An  alternative  picture,  that  a  target-gas  interchange 
instability  causes  the  turbulence,  is  less  likely.  If  (his  was  the  case,  we  would  intuitively 
expect  the  background  gas  to  slow  the  target.  Since  this  does  not  happen  we  believe  the 
first  picture  above  to  be  more  probable:  But  a  definitive  answer  awaits  more  experiments 
and  calculation. 
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FIGURE  4.  Sample  bright-field  and  phase-contrast  time  sequence  of  the 

turbulent  region.  The  target  is  20/im  CH  and  the  background  gas  is 
5  torr  of  nitrogen.  Laser  is  incident  from  the  bottom. 


FIGURE  3.  Trajectory  of  turbulent  region  for  die  case  of  a  9  fan  thick  CH  target 
irradiated  with  5  terawatts/cm1  in  5  ton  of  nitrogen.  The  turbulent 
region  moves  at  Mach  70. 


Nevertheless,  die  background  gas  does  mix  with  the  target  "piston",  and  is, 
therefore,  also  turbulent.  This  can  be  seen  by  examining  photographs  of  visible  light 
emission  from  the  turbulent  region  and  comparing  them  to  shadowgraphy  pictures  taken  at 
the  same  time  on  the  same  shot.  In  front  of  the  camera  is  a  narrow  band-pass  (50A 
FWHM)  interference  filter  which  passes  the  5001 A  Nu  line  from  the  background  nitrogen 
gas.  In  general,  the  overall  shape  of  the  emission  region  matches  closely  with  the  turbulent 
region  seen  by  shadowgraphy.  But,  small  scale  HMOOpm  structures  seen  in  the 
shadowgraphs  are  washed  out  in  the  emission  pictures,  probably  by  integration  of  die 
emission  along  the  line-of-sight  through  the  turbulent  volume:  The  smallest  emission 
structures  have  a  size  of  about  20Qpm.  (Camera  resolution  is  50fm  in  the  target  plane.) 
Limb  brightening  causes  the  outer  edges  of  the  emitting  region  to  be  brighter  than  the  inner 
regions.  Figure  6a  shows  a  sample  emission  picture  and  figure  6b  a  time  integrated 
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spectrum  from  the  turbulent  region  No  spectral  lines  from  highly  ionized  states  of  are 
observed.  From  this  we  estunate  that  the  temperature  in  the  turbulent  region  is  about  10  eV 
or  less. 


To  further  test  the  relationship  between  turbulence  and  background  gas  we  varied 
the  gas  pressure  (0.03  to  10,000  mTorT)  and  the  gas  species  (He,  N,  Xe).  Interestingly, 
below  200  mTorr  of  nitrogen  the  turbulence  could  not  be  seen  in  our  emission  and 
shadowgraph  diagnostics.  Turbulence  was  never  observed  in  high  vacuum.  This  may 
mean  that  the  shadowgraphy  diagnostics  are  sensitive  to  features  of  the  background  gas  • 
perhaps  density  gradients  in  the  swept  up  gas  -  and  that  collision  with  the  background  gas 
raises  the  temperature  so  that  the  region  emits  enough  to  be  seen.  This  would  have  to  be  the 
case  if  the  hypothesis  that  the  target  is  turbulent  without  the  presence  of  gas  is  correct 
When  the  turbulence  is  seen  it  does  not  seem  to  be  affected  by  the  type  or  the  pressure  of 
gas.  However,  detailed  spectral  information  has  not  yet  been  quantified. 


FIGURE  6.  (a)  Sample  emission  photograph  and  (b)  sample  emission  spectrum 

of  the  turbulent  region. 
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Power  Smctral  Iharjili-risliis  Nonuniformities  in  the  turbulent  region 

t-jn  be  quantified  hv  measuring  their  spatial  frequency  spectrum  Although  such  a 
spectnim  cannot  uniquely  identify  a  turbulence  mechanism,  any  theory  that  does  not 
reproduce  the  spectnim  can  be  eliminated.  In  addition,  such  spectra  are  used  by  designers 
of  systems  that  must  operate  in  turbulent  surroundings.  The  challenge  for  experiments  is  to 
relate  the  observed  structures  to  some  underlying  physical  quantity  such  as  density,  energy, 
index  of  refraction,  etc.  in  our  work  this  is  being  done  using  phase-contrast  shadowgraphy 
and  holographic  interferometry". 


Phase-contrast  shadowgraphy  is  an  application  of  the  Zemikc  phase-contrast 
microscope"  to  plasma  research.  The  method  is  described  with  modem  Fourier-transfotm 
mathematics  by  Cutrona  et  al  *  Presby  and  Finkelstein11,  we  recently  discovered,  were  the 
first  to  use  it  in  plasma  research.  The  essential  aspects  of  the  method  may  be  easily 
understood  by  considering  the  following  simplified  argument:  Imagine  the  laser  probe 
(Fig.  3)  to  be  a  planar  electric  wave  of  amplitude  E.  Assume  that  turbulence  alters  the 
phase,  but  not  the  amplitude,  of  the  probe  so  that  after  the  probe  traverses  the  turbulence  its 
electric  field  is  given  by  t 

(1)  E  ■  E,exp  { 1(9  +  fl(x,y) } 

■  E,exp(ip)  |  1  +  X(i0)*/n!  |. 

Here  9  is  an  averaged  spatially  independent  part  of  the  phase  shift  and  0(x,y)  is  a  spatially 
varying  part  (which  need  not  be  less  than  1).  The  probe  is  imaged  onto  a  camera  with 
dual-lens  astronomical  telescope  as  shown  in  Fig  7a.  The  first  lens  of  this  telescope 
Fourier  transforms  the  probe's  electric  field  at  its  focal  point.  The  second  lens  picks  up  the 
Fourier  transformed  electric  field  and  Fourier  transforms  it  again  at  the  location  of  die 
camera.  The  camera  photographs  the  square  of  the  twice  Fourier  transformed  electric  field. 
Mathematically,  the  first  Fourier  transform  is  given  by 

(2)  f[E]  -  E,  exp(ip)  |  F[I]  +  FT  I  (i$)*M!] ) ; 
and  the  second  Fourier  transform  by 

(3)  F[FTEJ  1»  Eo  exp(fp)  (  FIFTH  J+FTFTl  (i*)'/h!J] ) ; 
and  the  image  on  the  film  by 

(4)  Imaged)  =  |  FIFIE]]}1 

(Magnification,  time,  and  axial  direction  in  the  phase  terms  are,  for  simplicity,  not  included 
here.) 

Now,  the  FT1]  term  in  equation  2  is  the  transform  of  the  undisturbed  pan  of  the 
probe.  Physically,  h  is  a  spot  of  light  on  the  lens  axis.  This  spot  can  be  manipulated, 
thereby  modifying  the  photograph  given  by  equation  4. 

For  example,  without  any  manipulation  the  photograph  is  given  by 

(5)  lmage(x>y)  =  E^E,' . 

This  is  the  so  called  bright-field  shadowgraph.  In  the  absence  of  absorption  it  is  simply  the 
image  of  the  original  probe.  With  absorption,  h  forms  the  shadow  of  the  absorbing  regie 
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If  the  central  spot  is  eliminated  the  photograph  is  given  by 

(6)  Image(x.y)  *  EJE,,'  |  8\-x,-y)  +  ...|, 

which  is  a  Schlieren  shadowgraph  sensitive  to  the  square  of  phase  fluctuations. 

If  the  central  spot  is  shifted  by  a  quarter  wavelength  the  photograph  is  given  by 

(7)  Image(x,y)  ■  EeE,’  (  1  +  20(-x,-y)  4- ...  J  , 

which  is  a  phase -contrast  shadowgraph  sensitive  to  the  first  power  of  phase  fluctuations. 


From  phase  fluctuations  we  determine  fluctuations  in  the  index  of  refraction 
N(x,y,z)  integrated  along  the  axial  direction  using 

(8)  /N(x.y^)  dz  ■  A^j0(x,y)/2n 

where  is  the  probe  wavelength.  Assuming  free  electrons  of  density  n.(x,y,z) 
dominate  the  index  of  refraction  we  get 

a  i  -t 

(9)  JWxom)  dz  »  4>(x,y)  (A^e  /m,c )  . 


The  phase-contrast  method  is  subject  to  the  following  caveats: 


•  Significant  absorption  or  probe  nonuniformity  will  invalidate  the  commutation 
in  equation  2. 

•  As  a  practical  matter,  one  cannot  eliminate  FTl]  by  itself  One  also  eliminates 
the  DC  and  low  frequency  components  in  the  terms  containing  6.  Therefore, 
measurements  of  low  frequency  fluctuations  are  not  accurate  . 

•  If  phase  fluctuations  larger  than  n  occur  then  equation  7  cannot  be  uniquely 
inverted.  This  causes  aliasing  and  spectrum  distortions.  Therefore,  any  phase 
excursion  greater  than  it  should  not  occur  often  snd/or  be  at  wavelengths  not  of 
interest  to  the  analysis. 

•  Optical  distortions  are  neglected 


One  nice  feature  of  the  method  is  that  the  probe  is  imaged  right  after  leaving  the 
turbulent  region.  Consequently,  Fresnel  diffraction  effects,  which  would  occur  if  the  probe 
were  to  propagate  a  long  distance,  do  not  trouble  us  here.  Fresnel  oscillations,  which 
distort  phase  measurements  in  atmospheric  experiments13  are  not  a  problem. 


Fig.  7b  shows  how  the  method  is  implemented  in  practice.  The  probe  beam  is  split 
into  two  pans  by  a  Wollaston  prism  after  leaving  the  turbulent  region.  One  part  is  imaged 
as  a  bright-fleld  shadowgraph,  while  the  other  pan  is  passed  through  a  phase  filter 
producing  the  phase-contrast  shadowgraph  of  equation  7.  The  phase  contrast  filter  is  a  flat 
(A/10)  fused  silica  window  with  a  1mm  diameter  by  2859  angstrom  thick  (1/4  A)  mesa  on 
its  surface.  The  film  is  calibrated  by  pie-exposing  a  pan  of  it  to  the  probe  light  through  a 
Kodak  step  wedge. 
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FIGURE  7.  (a)  Schematic  of  the  phase-contrast  concept,  (b)  Implementation  of 

the  phase-contrast  method. 


Simultaneous  bright-field  and  phase-contrast  shadowgraphs  verify  that  absorption 
or  probe  nonuniformity  effects  ate  not  important  in  the  region  of  interest.  Also, 
photographs  taken  without  turbulence  or  plasma  being  present  let  us  balance  the  two  split 
beams  of  the  probe,  and  verify  that  extraneous  effects  from  the  mask  or  the  optics  are 
minimal.  The  diameter  of  the  mesa  inqilies  that  spatial  wavelengths  longer  than  300pm 
(«.  Ay-JLJlr _ )  are  not  measured  accurately. 


Figure  8  shows  typical  bright-field  and  phase-contrast  shadowgraphs.  The 
bright-field  image  allows  a  ahadow  of  the  uninadiated  part  of  the  foil,  and  a  shadow  of 
slow  material  ejected  by  shocks  and  heat  leakage  from  the  laser  spot.  Such  slow  material 
plays  no  role  in  turbulence  formation.  Also  visible  is  a  segment  of  the  ablation-side  blast 
wave  which  wrapped  itself  around  the  rear  of  die  foil.  This  blastwave  scatters  light  out  of 
the  optical  path  so  that,  for  practical  purposes,  it  behaves  as  an  absorbing  object.  The 
turbulent  region  is  barely  visible  in  a  few  tiny  spots  in  the  bright-field  image  but  is  very 
visible  in  the  phase-contrast  shadow,  proving  thru  the  turbulence  is  primarily  a  phase 
object.  The  phase-contrast  shadow  is  also  rich  in  other  phase  disturbing  features  such  as  a 
weak  shock  reflected  from  the  foil  holder,  and  sound  waves  outside  die  boundaries  of  the 
turbulence. 


The  phase -contrast  photographs  are  digitized  and  transformed  into  equivalent  phase 
images  on  a  computer  using  the  mathematics  outlined  above.  Then  die  spatial  power 
spectral  density  is  determined  from  the  expression 

(10)  Powerfic.) «  X,[Abs(jT(x)  0(x,y) exp(-ik,x) dx)  ] , 
where  T(x)  is  an  80%  Tulcey  window**. 
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FIGURE  8.  A  bright-fteld  and  a  phase -contrast  shadowgraph  of  a  10pm  A1  foil 
irradiated  with  4.5  twans/cm:  in  5  torr  of  nitrogen.  The  pictures 
were  taken  2 19ns  after  the  peak  of  the  laser  pulse.  Laser  is  incident 
from  the  left. 


A  typical  one-dimensional,  spatial  power  spectrum  is  shown  in  Fig.  9.  Most 
turbulent  power  spectra  can  be  characterized  by  three  quantities:  an  'outer  tcalelength'  at 
which  energy  is  pumped  into  the  system,  the  'inertial  range'  of  wavelengths  at  which  energy 
cascades  from  huger  to  smaller  scales,  and  an  ‘inner  scalelength'  at  which  energy 
out  of  die  system.  In  our  case  the  outer  scalelength  is  about  1mm  which  is  similar  to  the  \ 
diameter  of  die  focal  spot  (However,  this  number  is  not  very  reliable  as  explained  above.) 
The  inertial  range  has  a  k.4,1  scaling.  An  inner  scalelength  is  not  observed  at  20pm  or 
larger  wavelengths.  Inner  scalelengths  smaller  than  20pm  are  not  unreasonable  for  the 
estimated  Reynold’s  number  in  this  How.  The  scaling  of  the  inertial  range  is  independent 
of  scan  direction,  reflecting  symmetry  in  the  fine  turbulent  structure. 


A  two-dimensional,  k-space  power  spectrum,  PQk,  Jc,),  can  also  be  calculated.  We 
find  that  Pdc,Jc,«0) -k.^and  P(k,«0JcT)  =  k, "  P(k>«  /(Pflsjgk  d*.  where  *  is  the 
azimuthal  angle  in  k^k,  space,  varies  as  k*'.  The  various  scalings  are  all  consistent  for 
spatially  symmetric  turbulence. 


We  have  yet  to  show  that  unacceptable  aliasing  is  not  a  problem  with  the 
phase-contrast  measurements.  This  can  be  tested  with  holographic  interferometry  which 
measures  absolute  phase  disturbances  caused  by  phase  objects,  and  can  also  provide  an 
independent  determination  of  the  power  spectrum.  Such  intetferograms  have  been  taken 
with  a  variety  of  background  fringe  spacings.  Interferograms  with  fringe  spacings  of  1mm 
and  360pm  show  the  outlines  of  the  turbulent  region  dearly.  However,  the  fringes  cannot 
be  followed  into  the  turbulent  region  itself  since  they  have  a  coarser  scale  than  the 
turbulence  fluctuations.  But  fringes  in  interferograms  with  70pm  spacing  (Fig.10)  can  be 
followed  through  most  of  the  turbulence.  A  visual  analysis  shows  that  aliasing  is  probably 
not  significant,  but  a  quantitative  estimate  has  yet  to  be  made.  The  mathematics  for 
extracting  two-dimensional  phase  information  from  holographic  interferograms  of  laminar 
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shot. 


FIGURE  10.  Enlarged  section  of  a  holographic  interferogram  of  the  turbulent 
region.  Reference  fringes  are  spaced  7(tyrm  apart. 

phase  objects  have  been  discussed  by  Takeda25  and  Nugent1*.  We  are  currently  evaluating 
such  methods  for  accuracy  in  a  turbulent  case. 

Turbulence  control  Power  spectra  alone  are  insufficient  to  unambiguously 

identify  the  responsible  turbulence  mechanism.  Indeed,  we  have  not  proven  that  the 
nonunifonnities  are  turbulent  in  the  sense  that  energy  is  cascading  from  larger  to  smaller 
scales,  or  whether  they  are  a  frozen  eddy  structure  with  a  spectrum  determined  by  a  linear 
instability.  By  manipulating  the  turbulence  we  may  get  further  clues  into  its  nature.  For 
example,  we  can  try  to  turn  the  turbulence  on  and  off,  or  to  control  its  power  spectrum.  If 
the  Rayleigh-Taylor  instability  in  the  target  is  responsible  for  the  turbulence,  perhaps  we 
can  accelerate  a  stable  target  and  see  if  the  turbulence  disappears.  We  have  accomplished 
some  of  these  goals  in  a  few  different  ways.  For  example: 
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Burn  Through:  By  focusing  the  laser  beam  to  a  tight  spot  we  can  raise 

its  irradiance  so  that  the  entire  target  is  heated  by  x  rays,  fast  electrons  and 
explodes.  There  is  no  ablative  acceleration  and  the  Rayleigh-Taylor  instability 
is  stabilized.  When  we  irradiated  a  target  under  such  conditions  turbulence  did 
not  occur.  Instead,  we  observed  an  expanding,  hemispherical  front  similar  to 
the  blast  wave  normally  associated  with  the  target's  laser-side  (Fig.  11).  Long 
wavelength  nonuniformities.  similar  to  weak  regions  on  the  surface  of  an 
expanding  balloon,  replaced  the  fine-scale  turbulent  structure  associated  with 
ablatively  accelerated  targets.  This  observation  supports  the  turbulent  target 
model. 


ion 
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FIGURE  11.  Schlieren  shadowgraphs  of  a  9pm  CH  foil  in  5  ton  of  nitrogen 

irradiated  by  a  tightly  focused  laser  beam  to  320  terawatts/cm1.  Note 
lack  of  turbulence.  Laser  is  incident  from  the  left. 


Structured  targets:  The  Rayleigh-Taylor  instability  in  flat  targets 
develops  from  natural  material  imperfections  or  from  nonuniformities  in  the 
laser  beam.  The  wavelengths  which  grow  can  be  controlled  by  purposefully 
imposing  a  large  mass  perturbation  on  the  target  -  for  example,  by  grooving  its 
surface*’  *.  Now,  if  the  tutbulent  spectrum  is  determined  by  the  linear 
Rayleigh-Taylor  instability  then  we  should  be  able  to  alter  the  turbulent  power 
spectrum  by  changing  the  wavelengths  at  which  Rayleigh-Taylor  grows.  To  test 
this  we  accelerated  targets  perturbed  with  a  lOQpm  wavelength  groove  and 
compared  the  resulting  turbulent  power  spectra  to  those  of  non-perturbed 
targets.  We  found  that  at  early  (<10ns)  times  a  100pm  structure  was  visible  in 
the  shadowgraphs.  But  fully  developed  turbulence  at  200ns  had  the  same 
spectrum  as  that  for  flat  targets.  This  supports  the  argument  that  we  are 
observing  turbulence  and  not  a  frozen  eddy  structure. 
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SUMMARY 


This  paper  describes  our  studies  of  turbulence  produced  in  a  background  gas  by 
objects  accelerated  to  very  high  velocities.  We  have  developed  methods  to  quantify 
t  features  of  the  instability,  in  particular  the  power  spectrum  of  the  electron  density.  Also, 

we  have  measured  the  behavior  of  the  turbulence  under  many  different  circumstances  and 
determined  its  power  spectrum.  With  these  methods  we  are  studying  the  nonlinear 
hydrodynamics  of  a  system  which  duplicates  in  the  laboratory  many  astrophysics] 
properties.  We  intend  to  extend  our  methods  to  the  study  of  hydrodynamics  at  less  exotic 
conditions  but  without  the  complications  of  more  classical  experimental  methods. 
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